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I.  INTRODUCTION 


A.  Background 

The  calculation  of  radiant  intensity  from  rocket  and  jet-aircraft 
exhaust  plumes  requires  a priori  knowledge  of  the  pressure  p, 
temperature  T,  and  active  molecular  species  concentration  c distri- 
butions (pTc  distributions)  throughout  the  whole  plume  flow  field.  Direct 
measurements  of  these  distributions,  particularly  near  the  engine  nozzle 
exit  plane,  are  often  impractical  due  to  the  harsh  environment  of  the  plume. 
In  addition,  direct  measurements  are  intrusive  and  often  disturb  the 
flow  field  to  such  an  extent  that  meaningless  re  suit  are  obtained. 
Accordingly,  most  plvime  flow-field  determinations  are  based  on  computer 
modeling. ' ' Computer  codes  that  generate  the  downstream  properties 
of  the  plume  often  require  the  radial  pTc  profiles  at  the  nozzle  as  initial 
boxindary  conditions  for  the  solution.  These  profiles  are  calculated  with 
still  other  codes  from  the  thermodynamics  of  the  propellant  combustion, 
nozzle  expansion,  and  other  engine  characteristics.  Although  many  of  these 
codes  are  quite  sophisticated  and  based  on  sound  physical  and  chemical 
principles,  there  is  a decided  lack  of  experimental  verification  of  their 
predictions. 

In  order  to  assess  the  accuracy  of  those  code  predictions,  a series 
of  experimental  programs  to  infer  radial  pTc  profiles  from  nonintrusive 
measurements  of  plume  emission  and  absorption  (E/A)  characteristics  was 
initiated  by  the  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL).  All  of 
the  daca  for  these  programs  were  obtained  from  engine  test  firings  at  the 
Arnold  Engineering  Development  Center  (AEDC)  facility  by  ARO,  Inc. 

Some  results  and  analyses  of  data  from  the  first  two  programs  [High  Alti- 
tude Plume  Radiation  Program  (HAPRAP)  and  Emitted  Radiation  from 
Special  Engines  (ERASE)  program]  are  reported  in  Ref.  3 and 
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Refs.  4 and  5,  respectively.  The  present  work  was  performed  for 
the  current  Nozzle  Exit  Radiation  Determination  (NERD)  program*  and  con- 
cerns the  deduction  of  plume  pTc  properties  from  measured  E/A  data  by 
mathematical  inversion  of  the  equations  of  radiative  transfer. 

B.  Review  of  Emission/Absorption  Inversion 

The  basic  principle  of  emission/absorption  inversion  as  applied  to 
gases  can  be  understood  by  considering  a \3niform  optical  path  of  length 
L,  pressure  p,  mole  fraction  of  active  gas  c,  and  temperature  T.  For 
given  p and  L,  the  two  \inknowns  are  c and  T,  and,  consequently,  a 
measurement  of  at  least  two  independent  properties  of  the  gas  sample 
will  be  required  to  determine  these  values.  The  radiance  and  absorptance 
(or  transmittance)  of  the  gas  sample  at  a given  wavelength  (or  within  a 
small  wavelength  interval  centered  about  the  wavelength)  are  such 
independent  properties.  The  temperature  can  be  determined  from  the 
measured  radiance  and  absorptance  by  application  of  Kirchhoff's  law, 
and,  once  the  temperature  is  known,  the  value  of  c can  be  derived  from 
the  absorptance  measurement  by  use  of  appropriate  transmission  formulas. 

The  same  principle  can  be  applied  to  nonuniform  optical  paths,  but 
its  application  is  seriously  compounded  by  the  existence  of  both  tempera- 
ture and  species  concentration  gradients.  The  general  method  of  solution 
involves  the  assumption  that  the  nontmiform  path  can  be  divided  into  a 
number  of  intervals,  and  that  the  properties  within  each  interval  can  be 
considered  as  constants.  If  n intervals  are  used,  the  number  of  unknown 
quantities  is  now  2n  (T  and  c of  each  interval).  Thus,  Zn  independent 
quantities  must  be  measured  in  order  to  extract  a solution. 

The  two  main  approaches  to  measuring  the  requisite  number  of 
independent  quantities  are  to  measure  the  radiance  and  absorptance  at 
either  n different  spectral  positions  (multi-color),  or  for  n independent 
lines  of  sight  (multi-position)  through  the  gas  sample.  For  an  axisym- 

*W.  Bertrand,  T.  McRay  and  E.  Kiech,  Nozzle  Exit  Plane  Radiation  Diag- 
nostics of  Three  Transtage  Engine  Injector  Configurations,  AEDC  Tech- 
nical Report,  to  be  published  (1978). 
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metric  radiation  source  such  as  a rocket  plume  near  the  exit  plane,  either 
of  these  approaches  is,  in  principle,  applicable.  Multi-color  measure- 
ments could  be  made  through  the  diameter  of  the  plume  or  multi-path 
measurements  at  a fixed  spectral  position  could  be  made  by  traversing 
across  the  lateral  extent  of  the  plume.  The  viewing  direction  for  either 
of  these  measurements  would  be  perpendicular  to  the  plvime  axis  so  that 
the  results  pertain  to  a single  axial  station  of  the  plume. 

Techniques  for  performing  both  multi-color  and  multi-position 
measurements  and  inversion  have  been  considered  extensively  with  regard 
to  the  problem  of  atmospheric  soianding  by  satellite -borne  sensors. 

Such  passive  techniques  are  able  to  yield  altitude  profiles  for  temperature 

and  concentration  of  H^O  and  O^.  Krakov,  Simmons  et  al. , 

(9)  ^ ^ (10) 

Buchele, ' ' and  Cutting  and  McC.  Stewart'  ’ discuss  mult-.-color  techni- 
ques as  applied  to  high  temperature  combusilcn  gas  pioblems.  The 
technique  of  Krakov  requires  that  radiance  measurements  be  performed 
at  spectral  positions  where  the  sample  is  described  by  strong  absorption 
since  his  band  model  radiation  formulation  for  transmittance  is  limited 
to  this  condition.  The  inversion  algorithm  employs  the  techniques  of 
successive  approximations  and  Lagrange  undetermined  multipliers.  The 
use  of  the  latter  technique  allows  a solution  to  be  obtained  for  the  over- 
determined system  of  equations  resulting  from  the  use  of  measurements 
at  more  than  n different  spectral  positions.  The  instabilities  that  often 
result  from  overdetermined  systems  and  measurement  errors  in  the 
input  data  is  balanced  by  choosing  m\iltiplier  values  which  stifficiently 
"smooth'' the  data  to  yield  convergent  iterations. 

The  techniques  of  Simmons  et  al.  and  Buchele  are  based  on  general 
band  model  formulations  that  do  not  require  the  limitation  to  strong 
absorption.  The  "inversion"  involves  the  assumption  of  a priori  species 
concentration  and  temperature  profile  functions  which  contain  a few 
adjustable  parameters.  The  parameters  are  varied  until  the  spectrum 
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generated  from  these  profiles  using  the  band  model  formulation  best 
matches  the  observed  spectrum.  Justification  for  the  use  of  a priori 
profile  functions  is  based  on  the  assumption  that  the  physical  restraints 
of  most  systems  will  indicate  a good  choice  of  profile.  A principal 
advantage  of  this  technique  is  that  the  problem  of  ntimerical  instabilities 
inherent  in  other  inversion  techniques  is  essentially  obviated.  Moreover, 
if  the  functional  form  of  the  profiles  is  chosen  general  enough  (i.  e. , it 
contains  many  parameters),  this  procedure  approaches  an  "inversion"  in 
the  usually  interpreted  sense  of  the  word. 

The  method  of  Cutting  and  McC.  Stewart  is  very  nearly  a direct 
application  of  atmospheric  sounding  techniques  to  high -temperature 
furnace  conditions. 

The  multi -position  technique  applied  to  high  temperature  combustion 
gas  systems  is  discussed  by  Brewer  and  Limbaugh.  ^ ^ ^ ^ In  this  approach, 
the  spatial  radiance  and  absorptance  cf  an  axisymmetric  cylindrical 
source  are  measured  for  a single  wavelength  interval  A X as  ?,  function 
of  transverse  distance  across  the  source.  The  measured  profiles  thus 
represent  results  for  all  the  possible  independent  lines  of  sight  through 
the  sample  (for  a perpendicular  viewing  geometry).  The  inversion  techni- 
que assumes  that  the  axisymmetric  source  is  composed  of  n homogeneous 
annular  zones.  The  edge  of  the  measured  profiles  thus  correspond  to  a 
line  of  sight  passing  through  the  single,  outermost  homogeneous  zone. 

The  center  of  the  profiles  corresponds  to  a line  of  sight  through  a series 
of  2n-l  homogeneous  segments.  By  starting  at  the  edge  of  the  measured 
profiles  and  working  toward  the  center,  an  inversion  can  be  effected  by 
iteratively  determining  the  T and  c of  each  new  zone  that  is  added  to  the 
lines  of  sight.  This  inversion  method  constitutes  the  primary  diagnostic 
tool  that  has  been  used  in  the  analysis  of  the  HAPRAP  and  ERASE  data 
and  is  presently  being  used  by  ARO,  Inc.  in  the  NERD  program. 
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C.  Scope  of  Present  Work 


The  present  work  is  an  in-depth  assessment  of  the  multi -position 
inversion  method  for  axisymmetric  cylindrical  pl\ime  sources.  The 
fundamental  goal  is  the  inversion  of  profiles  of  plume  radiance  and 
absorplance  obtained  in  transverse  scans  across  the  lateral  extent  of  test 
engine  plumes  near  the  nozzle  exit  plane  in  order  to  obtain  the  radial 
profiles  of  temperature  and  species  concentration,  and,  to  estimate  the 
accuracy  with  which  these  profiles  can  be  obtained.  The  principal 
radiating  species  and  spectral  regions  of  interest  are  H^O  (2.44-  2.  63//m) 
and  CO^  (4.  17  - 4.  55  //m). 

An  important  aspect  in  the  formulation  of  inversion  diagnostic  tools 
is  the  radiation  transport  model  employed.  The  present  work  uses  band 
model  methods  in  which  high  spectral  resolution  is  sacrificed  in  order  to 
gain  efficient  computation  of  mean  E/A  values  in  spectral  bandpasses 
encompassing  several  tens  or  hundreds  of  individual  vibration-rotation 
lines  of  a molecular  band.  Use  of  such  models  is  entirely  consistent 
with  the  spectral  resolution  (*~  20 -300  cm”  ) with  v/hich  experimental  E/A 
profiles  are  obtained  in  the  AEDC  measurements  programs.  The  band 
model  formulation  used  in  the  present  work  is  a state-of-the-art  model 
for  application  to  combustion  gas  problems.  The  three  principal  features 
that  distinguish  this  formulation  from  many  previously  used  formulations 
are:  (i)  account  of  Doppler  broadening  effects  in  low  pressure  plumes 
through  the  use  of  Yoigt  line  shape  models;  (2)  the  treatment  of  nontiniform 
source  conditions  (inhomogeneities  and  nonisothermalities)  with  approxi- 
mations better  suited  than  the  traditional  Curtis -Godson  (CG)  approximation 
to  the  sharp  gradients  encountered  in  many  plumes,  and  (3)  the  use  of  band 
model  parameters  that  are  consistent  over  the  temperature  range  from 
200  K to  3000  K.  A general  outline  of  the  radiation  model  is  presented 
in  Section  II  and  considerations  on  band  model  parameters  for  the  current 
applications  are  given  in  Section  V. 
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Once  the  equations  relating  the  E/A  profiles  to  the  radial  pTc  profiles 

are  defined  through  the  band  model  formulation,  a means  of  inverting  these 

radiation  transport  equations  must  be  constructed.  When  the  plume  source 

is  optically  thin,  the  transport  eqviations  reduce  to  Abel  integrals  and 

l\2.) 

standard  techniques  of  inversion  may  be  applied. ' ' For  nonthin  sources, 

iteration  methods  must  be  used.  For  the  present  work,  an  iterative  Abel 
inversion  algorithm  was  formulated.  The  procedure  is  initiated  with  an 
Abel  inversion  (regardless  of  the  s /urce  optical  thickness)  and  the  radial 
profile  results  iteratively  corrected  by  performing  Abel  inversions  on  the 
difference  profiles  obtained  by  subtracting  the  input  E/A  profiles  from 
intermediately  computed  E/A  profiles. 

A necessary  procedure  in  any  inversion  of  real  data  is  some  form 
of  data  smoothing.  Regardless  of  the  inversion  method,  experimental 
random  noise  on  the  input  E/A  profiles  is  amplified  in  the  inversion  process 
and  can  yield  radial  profiles  so  noisy  that  they  are  meaningless.  The 
smoothing  procedure  adopted  for  this  work  is  a straightforward  pre- 
smoothing of  the  data.  Both  the  iterative  Abel  inversion  algorithm  and 
the  presmoothing  method  are  presented  in  detail  in  Section  in. 

During  the  course  of  this  work,  a second  inversion  method  (the 
Phillips -Twomey  method)  that  provides  implicit  smoothing  of  the  input 
data  was  formulated  and  compared  with  the  iterative  Abel/pre  smoothing 
method.  The  method  was  not  found  to  be  significantly  better  than  the 
iterative  Abel  method.  The  formulation  and  one  comparison  are  presented 
in  Appendix  B of  this  report.  Aside  from  this  appendix,  all  of  the  work 
of  this  report  employed  the  iterative  Abel  inversion  scheme. 

This  report  is  divided  into  two  major  parts.  Part  1 contains  the 
formulation  of  the  general  radiative  (Section  II)  and  inversion  (Section  III) 
models  for  axisymmetric  plume  radiation  diagnostics.  The  incorporation 
of  these  models  into  a working  computer  code  (EMABIC)  is  described  in 
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Section  IV.  Part  2 of  the  report  is  an  assessment  and  application  of  the 
overall  formulation.  Band  model  parameters  appropriate  to  the  present 
applications  are  discussed  in  Section  V.  Examples  of  E/A  profile  pre- 
dictions from  given  pTc  profiles  for  a variety  of  engine  types  are  presented 
in  Section  VI.  The  inversion  characteristics  of  the  iterative  Abel  inversion 
method  are  assessed  in  Section  Vn,  and  the  effects  of  random  and  bias 
data  error  are  treated  in  Sections  Vin  and  IX,  respectively.  An  application 
of  the  formulation  to  the  inversion  of  actual  NERD  E/A  data  is  made  in 
Section  X.  A summary  and  conclusions  section  (Section  XI)  completes 
the  main  text  of  the  report.  Three  appendixes  treat  topics  that  do  not  fall 
easily  into  the  mainstream  of  the  report.  Considerations  on  the  selection 
of  optimum  spectral  bandpasses  for  future  E/A  measurements  programs 
is  presented  in  Appendix  A,  the  Phillips -Twomey  inversion  scheme  is 
described  in  Appendix  B,  and  a discussion  on  the  construction  of  band 
model  parameters  for  wide -band  application  is  given  in  Appendix  C. 
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PART  1 


band  model  radiation 
and  inversion  modei^ 


n.  BAND  MODEL  RADIATION  FORMULATION 


A.  Introduction 

Radiation  calculations  are  performed  entirely  within  the  statistical 
band  model  formulation  for  an  array  of  lines  distributed  in  strength 
according  to  the  exponential -tailed  inverse  distribution.  Derivations  and 
discussions  of  the  basic  models  and  approximations  that  account  for  non- 
uniformities along  optical  paths  and  Doppler  broadening  effects  are  given 
in  Refs.  13  through  16.  Only  the  equations  defining  the  models  and  methods 
of  calculation  are  given  here.  In  parts  B and  C of  this  section,  the  radia- 
tion formulations  for  a general,  nonuniform  optical  path  are  reviewed.  The 
coordinate  transformation  to  the  cylindrical  geometry  relevant  to  plume 
radiation  is  treated  in  part  D. 

B.  Radiance  and  Transmittance  Equations 

Consider  a general  optical  path  along  which  the  total  pressure  p(s), 
temperature  T(s)  and  concentration  (mole  fraction)  of  active  gas  c(s)  are 
known  as  a function  of  geometric  path  position  s.  The  total  length  of  the 
path  is  s^,  and  the  observation  position  is  s = 0.  The  fundamental  band 
model  para-  ie<ers  k(s)  (mean  absorption  coefficient),  6(s)  (mean  line 
spacing)  and  y(s)  (mean  line  width)  for  the  spectral  interval  Av  and  active 
gas  species  are  also  presumed  known  as  a function  of  s.  The  first  two 
parameters  are  functions  of  T(s)  only,  whereas  y(s)  may  depend  addition- 
ally on  c(s)  and  p(s).  Explicit  tabulations  of  these  band  model  parameters 
for  H^O  in  the  2.  5-pm  region  and  in  the  4.  4-nm  region  are  given  in 
Section  V. 

The  mean  value  of  radiance  in  the  interval  Av  about  v is  given  by 

®m  _ 

N = - V N*(8)^^ds  (1) 

0 
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where  N (s)  is  the  Planck  radiation  function  evaluated  at  T(s)  and  v.  t(s) 
is  the  mean  transmittance  in  Av  for  the  path  segment  between  s = 0 and 
general  position  s.  is  the  total  geometric  length  of  the  line  of  sight, 

and  no  source  exists  beyond  s . Within  the  statistical  band  model  formu- 
lation,  the  transmittance  is  given  by 

T(3)  = ( 

where  W(s)/6  measures  the  ratio  of  the  mean  equivalent  width  (total  inte- 
grated absorptance)  of  the  lines  in  Av  to  the  mean  line  spacing  in  Av.  The 
transmittance  derivative  is 


diM 

ds 


-T{s) 


1 dW(s) 
6 ds 


(3) 


For  all  lineshapes  and  nonuniformity  approximations,  the  equi /alent  width 
derivative  and  equivalent  • .’idth  are  determined  from 

= c{s)p(s)k(s)y(s)  (4) 


and 


s 

= ^c(s')p(s’)k(s’)y(s')ds'  (5) 

0 

where  only  the  derivative  function  y(s)  changes  with  model.  These  func- 
tions are  considered  in  Section  C.  With  the  relationships  (1)  through  (5), 
the  total  radiance  and  transmittance  of  the  optical  path  may  be  written 


N 


\ c(s)p{s)k(s)T{s)N  (s)y(s)ds 
0 


T 


exp 


s 

m 


- ^ c{s)p(s)k(s)y(s)ds 
0 


(6) 


(7) 
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In  order  to  simplify  the  notation,  define 


w(s)  = T{s)y(s) 

(8) 

K(s)  = c(s)p(s)k(s) 

(9) 

J(s)  = N*(8)K(8) 

(10) 

J(s)  and  K(s)  are  called  the  source  function  and  absorption  function. 

res- 

pectively.  In  terms  of  K, 

J and  w,  the  fundamental  relations  for  radiance 

and  transmittance  are 

s 

m 

N = ^ J(s)w(s)ds 

(11) 

0 

s 

m 

-in  ^ K(s)y(s)d8 

(12) 

C.  Derivative  Functions 

The  preceding  formulation  fully  defines  the  radiation  and  transmis- 
sion of  the  optical  path  once  the  function  y(s)  is  defined.  These  functions 
are  now  considered.  In  the  traditional  Curtis -Godson  (CG)  appro'iimation 
for  an  array  of  Lorentz  lines 


y(s)  = y[nx(s),  p(s)] 


(13) 


where 


with 


y(Ti.p)  = (2-p)^^+  (p-1)^ 


1(T1)  = y/  1 + 2ti‘-  1 
dl(Ti)  _ 1 


(14) 

(15) 

(16) 
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x(s)  = 


u(B)k  (s) 


Ols)  - 

- 3 Js) 

6 


u(s)  = \ c(s*)p(s*)ds* 


ke(s)  = ^ c(s')p(s')k(s')ds' 


3 (s)  = 1 ^ c(s')p(s')k(s')0(s')d£ 

u(s)k^(s)  - 

3(3)  . 2^ 


In  thi««  approximation,  the  relation  (5)  for  W{s)/6  can  be  integrated  in 
closed  form  to  give 

W(s)  , XT 

= — jp-Ilnx(s)J 

(15) 

In  the  more  sophisticated  derivative  (DR)  approximation  ' for 
Lorentz  lines. 


y(s)  = y[TTx(s),  p(s),  q(s)] 


where 


? r r ^ q dz 

y(T^.  P.q)  =-\  \ exp  2IL__  du  j 4.  ^2  (2! 

0 [ 0 I ^ P ^ 1 J 

A numerical  procedure  for  computing  this  function  is  discussed  in  Ref.  15. 
Briefly,  the  procedure  consists  of  rewriting  y('n,  p,  q)  in  an  alternate  form 
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y(i'l,  p,  q) 


(26) 


1 

= ~ yC'H.  p»  1)  + y(xz*^,  p,  i)dz 

0 

that  relates  y(T),  p,  q)  to  its  value  at  q = 1.  The  integral  term  of  this  alter- 
nate expression  is  evaluated  numerically  by  Gaussian  quadrature,  and  the 
function  7(^1,  p,  1)  is  approximated  by 


y 1 + 2ti’  I^p  + y 1 + 


(27) 


which  is  the  derivative  function  within  the  Lindquist-Simmons  (LS)  approx- 
imation. The  additional  parameter  q(s)  of  (24)  is  defined  as 


6 (s) 
q(s)  = •::: — 
6(6) 


(28) 


where 


s 

J (s)  ::  3 \ c(s')p(s')k(s’)6(B')ds’  (29) 

u(s)k^(s)  Q 

In  this  approximation,  W(s)/£>  must  be  computed  from  (5). 

Fcr  an  array  of  Doppler  lines,  the  formulation  analogous  to  the  CG 
approximation  is 

y(s)  = y[TTx(s),  p(s)]  (30) 

where 

y(Ti,  p)  = (2  - p)  + (p  - 1)  (31) 

and 


H(Ti)  = — ^ In 
n n 


1 + Tie 


-u 


du 


(32) 
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The  functions  H(rj)  and  dH('n)/dTi  are  obtained  by  interpolation  on  the  data 
of  Table  1,  from 


H(ti) 


-E 


irii 


n+1 


^ (n  + 1) 


3/?. 


n=0 


dH(n) 

dri 


>/  n + 1 
n=u  ^ 


for  small  r„  and  from  the  asymptotic  series 
H(n) 


dH(r,) 

dr; 


3'^ 

2 

r /~TT 


1 4.  , 7tt‘^/64Q  , 

1 ^ . T . . . 

(inri)  (inn) 


(fnT^) 


1/2 


tt^/24  _ 35^^^/1920  ^ 
(inri)^ 


(33) 


(34) 


(35) 


(36) 


for  large  n.  The  arguments  x(s)  and  p(s)  are  defined  by  (17)  through  (21) 
and  with  5(s)  defined  by 


e(s) 


(37) 


where  vis)  is  the  Doppler  line  width.  Within  this  approximation  for  the 
Doppler  line  shape  (which  will  be  referred  to  as  th«  CG  approximation  for 
identification),  W(s)/6  is  given  in  closed  form  by 


W(s) 

6 


H[Tix(s)] 


(38) 


The  derivative  approximation  for  the  Doppler  line  profile  is 


y(s)  = 


where 


(3b) 
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f. 

I- 

f 


I 


£ 


1' 


1 


»• 

i 

I 


Doppler  Line  Band  Model 


■n 

H(n) 

dH(ri)/dTi 

0.  1 

0.09665 

0.  93460 

0.  2 

0. 18722 

0. 87828 

0.  4 

0.  35325 

0. 78593 

0.  6 

0. 50289 

0.  71310 

0.8 

0.63940 

0.  65398 

1. 0 

0.76515 

0.  60490 

2.  Q 

1.28138 

0. 44564 

A.  0 

2. 00567 

0. 29978 

6.  0 

2. 52794 

0. 22935 

8.0 

2.94136 

0. 18717 

10 

3.28568 

0. 15883 

20 

4.47964 

0. 092598 

40 

5.84661 

0. 052199 

60 

6.71916 

0.  036911 

80 

7. 36866 

0. 028752 

100 

7.88910 

0. 023642 

200 

9.  59377 

0.  012759 

400 

11.4233 

0.  006813 

600 

12. 5478 

0.  004702 

800 

13. 3689 

0.  003609 

1000 

14. 0187 

0. 002937 
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y{Ti.  P)  = ^ 

'^■6[l  + Tie-P  ==  j 

y('n,  p)  is  the  derivative  function  for  Doppler  lines  in  the  LS  approxima- 
tion, and  the  modification  of  the  optical  depth  argument  from  itx(s)  to 

TTx(s)/q(s)  result?  from  an  application  of  the  mean-line  (ML)  approxima- 

( 1 5) 

tion  for  Lorentz  lines  to  Doppler  lines.  y(ri,  p)  is  obtained  by  interpo- 
lation on  the  data  of  Table  Z,  from 


y(Ti,  p)  = 


n=0  Y 1 + np 


for  small  q,  or  from 


y(r,  p) 


Y(1  + n)[l  - T^(p^  - 1)] 


for  p < 0.  5 (error  < 17«).  For  t'>  10  , y(r|,  p)  is  set  to  zero. 

For  an  array  of  Voigt  lines,  two  approximations  are  used.  The  first 
is  an  heuristic  adaptation  of  the  approximation  devised  by  Rodgers  and 
Williams  (RW)  for  isolated  lines.  The  derivative  function  is 

y(s)  = A(s)yj^(s)  + B(s)yp(s)  + C(s)  (4: 


where 


A(s)  = 


W^(s)/6 

W(s)/5 


\^W^s)/6/| 


B(z)  = - 


>d(8)/6 


W{8)/6 


Wj8)/6 

W^{3)/6 


Table  2 • Equivalent  Width  Derivation  function  for  Statistical  Do 


(46) 


C(5) 


W^,(s)/6  /Wj  (s)/6 
W(s)/6  \w^(s)/6 


1 

ri 

2 

" ___ 

2 

■ f 

W(s)  _ \ 

Wj_^(s) 

Wd(s) 

Wj^(s)  Wp(s)  /W^(s) 

f - V 

6 

L 6 J 

L 6 6 / 6 J 

and 

= u(s)k^(s)  (48) 

Wj^(s)/6  and  yL(s)  are,  respectively,  the  equivalent  width  function  and 
derivative  function  evaluated  as  though  the  lines  were  pure  Lorentz  lines 
(in  either  the  CC  or  DR  approximations).  Wjj(s)/6  and  y£)(s)  are  the 
respective  functions  evaluated  assuming  a pure  Doppler  profile.  This 
approximation  is  accurate  to  better  than  ~ 10%. 

During  the  course  of  this  work,  a second  approximation  was  con- 
structed. It  is  essentially  the  same  as  the  Voigt  modei  employed  in  the 
NASA  plume  radiation  code,  ' ' but  uses  the  present  methods  for  comput- 

ing the  curves  of  growth  W^(s)/6  and  Wj^(s)/6  rather  than  che  methods 
prescribed  in  the  NASA  code.  The  derivative  function  y(s)  is  given  by  (43) 
where  now 


A{s) 


B(s) 


W Js)/6 

1 - 

W^(s)/6 

2 

W(s)/6 

W^(s)/6_ 

Wj^(s)/6’ 

1 - 

Wj^{s)/6 

2 

W(s)/6 
• « 

W^(s)/6 

(49) 


(50) 
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W(8) 

6 


1//y(s)‘ 


1 - 

'w^(s)/6‘ 

2 

“4- 

4. 

1 - 

’Wj^(8)/6' 

2 

W^(s)/6 

X 

W^(8)/5 

. 

(51) 


(52) 

(53) 


This  approximation  is  accurate  to  better  than  ~ S%. 


D.  Transformation  to  Cylindrical  Coordinates 

Consider  the  geometric  configuration  of  Fig.  1 in  which  measure- 
ments of  radiance  and  absorptance  are  made  on  a cylindrical  radiation 
source  at  a distance  z from  the  cylinder  diameter.  The  radiation  formu- 
lation has  been  developed  so  far  in  terms  of  the  optical  path  variable  s. 
We  wish  to  transform  from  s to  the  radial  coordinate  r with  z retained  as 
a parameter.  This  transformation  is  achieved  through 


■ ~"V'f»5»wr -*■.?■■» 


s (z)  = 2 - z^ 


Application  of  the  radiation  band  model  requires  the  evaluation  of 
various  integrals  over  the  optical  path  coordinate  s.  It  is  these  integrals 
that  are  transformed  from  s to  r.  Let  g(z,  s)  be  a general  function  of  e 
for  the  iine-of- sight  chord  at  z,  and 


f(z,  s}  = ^ g(z,  s*)ds' 


Transformation  of  this  general  integral  form  to  cylindrical  coordinates 
with  (54)  through  (57)  gives 


F (z,  r)  10  s 0 

F(z,  r)  = . 

F (z,  r)  cp  s 0 


where 


F“(2,  r)  = \G'(z,r') 


/“,2  ? 
V r'  - z 


^ r /-O,  r'dr'  r ,v  ’•'dr* 

r;  = 2 \ G (z,  r') " \ G (z,r')— 

/ .2  2 J~T~~ 

z ''r'-zr  Vr'- 


,2  2‘ 
r'  - z 


F(z,  r)  = f{z,  8(r)j 

s (z)  . — •, yl 

G (z,r)  = g z,  — r - z 

r«+/  » . v 2 2 1 

G (z,  r)  = g z,  — ^ v r - z 

J 

G^(z,r)  = "I  |G“(z,r)  + G^(z,  r)j 


band  model  parameters  are  then  also  functions  or  r only,  k = k(r),  6 = 6(r), 
V = y(r)  and  0 = 0(r),  as  are  the  absorption  (9)  and  source  (10)  functions, 

K = K(r)  and  J = J(r),  respectively.  Since  the  integrals  appearing  in  (19) 
through  (21)  and  (29)  for  the  optical  depth  u and  path  averaged  band  model 
parameters  k^,  0^  and  6^,  respectively,  all  involve  integrands  which  are 
functions  of  r only,  the  restricted  transformation  equations  (67)  through  (69) 
can  be  applied  to  compute  the  "positive"  and  "negative"  components 

u^(z,  r)  u"(z,r) 

k^(z,  r)  k^{z,  r) 

0^(z,  r)  3g(z,r) 

r)  6~(z,  r) 

from  which  the  parameters 

•* 

X (z,  r)  X (z,  r) 

p''‘(z,r)  p‘(z,r) 

q’*’(z,r)  q'(z,r) 

can  be  found  from  (17),  (18)  and  (28),  respectively.  These  latter  para- 
meters define  two  components  for  the  derivative  function  y(x,  p,q) 


y 


(z,  r) 


^ db  ^ 

y{x  (z,  r),  p (z,  r),  q (z,  r)] 


and  the  average 


(71) 


y®(z,  r)  = I [y^(z,  r)  + y"(z,  r)] 


(72) 


With  y (z,  r)  determined,  the  equivalent  width  (5)  is  determined  by  appli- 
cation of  the  general  transformation  aquations  (59)  through  (61)  [since  the 
integrand  is  now  a function  of  both  r and  z through  the  appearance  of 
y*(z,  r)  in  the  integrand]  to  obtain  the  two  components  of  W/6 


W^(z,r) 

6 


(73) 
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These  components  are  used  with  (2)  to  get  the  two  components  of  trans- 
mittance 


III.  'NVEPSION  MODELS 


A,  General  Considerations 


The  inversion  procedures  considered  in  this  section  are  based  on 
the  radiative  transfer  equations  (78)  and  (79)  for  transverse  radiance  and 
transmittance,  respectively.  With  the  definition 


these  equations  are 


W(z)  = -fnT{z) 


R 


N(z)  = 2 ^ J(r)w^(z,  r) 

/“2  ? 

z ^ r - z 


W(z)  = 


R 

2 \ K(r)y°(z,  r) ^ 


/”2  2 

r - z 


^80) 


(81) 


(82) 


The  basic  aim  of  inversion  is  to  determine  the  radial  functions  J(r)  and 
K(r)  given  the  transverse  functions  N(z)  and  W(z),  When  this  has  baen 
accomplished,  the  determination  of  the  radial  profiles  of  c,  p and  T pro- 
ceeds directly  from  the  definitions  (9)  and  (10) 


K(r)  = 

c(r)p(r)k(r) 

(83) 

J(r)  = 

K(r)N*(r) 

(84) 

'T’he  ratio  J(r)/K(r)  is  the  Planck  function  N (r).  Thus,  the  radial  temper- 
ature profile  is  determined  as 


T(r) 


in  C 


K(r) 
1 J(r) 


(85) 


where  Cj  = 1. 129  x lO”^^  cm^  W/cm^-sr-pm  and  = 1.4389  v cm  K. 
With  T(r)  known,  k(r)  is  known.  Then,  tlie  partial  pressure  c(r)p(r)  is 
obtained  from 
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c(r)p{r)  = 

k(r) 


(86) 


Note  that  the  inversion  solves  only  for  the  product  c(r)p(r)  and  not  the  two 
profiles  c(r)  and  p(r)  individually.  Even  though  the  product  is  usually  the 
final  result  desired  in  the  inversion,  it  is  significant  that  the  individual 
profiles  cannot  be  determined.  In  all  but  the  Abel  inversion,  an  iterative 
procedure  is  employed  wherein  (81)  and  (82)  are  alternately  used  for 
inversion  and  for  generating  intermediate  transverse  profiles.  While  the 
inversion  gives  only  c(r)p(r),  the  individual  variables  are  required  in  the 
profile  generation  through  the  calculation  of  y (see  Section  V).  The  math- 
ematical reason  for  this  problem.  Is  that  by  inverting  only  two  independent 
profiles,  that  is  N(2)  and  W(2),  we  are  allowed  to  obtain  only  two  radial 
profiles,  that  is  T(r)  and  c(r)p(r).  In  practice,  inversions  are  carried 
out  with  a pressure  profile  which  is  assumed  to  be  close  to  the  true  pres- 
sure profile.  A test  inversion  is  presented  in  Section  IX,  in  which  p(r)  is 
varied  about  a known  true  profile.  The  inversion  results  for  cvr),  of 
course,  vary  with  p(r),  but  the  partial  pressure  profile  c(r)p(r)  is  rela- 
tively insensitive  to  the  variation. 


B.  Abel  Inversion 


For  a source  that  is  optically  thin,  the  derivative  function  y(s), 
the  transmittance  t(s),  and  the  product  w(s)  = y(s)T(s)  approach  unity.  In 
this  limit,  the  radiative  transfer  functions  (81)  and  (82)  approach 


N(2) 


R 

2 ^ J(r) 


rdr 


(87) 


W(z) 


rdr 


In  both  cases,  the  equations  are  of  the  Abel  -.ntegral  form 


(88) 
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rdr 


(89) 


f(z)  = 


R 


for  which  the  well-known  inversion  is 


g(r) 


1 dF(r) 
2 nr  dr 


■i 


(90) 


F(r)  = 


(91) 


Thus,  solutions  for  J(r)  and  K(r)  are  iminediateiy  obtained  by  application 
of  (90)  and  (91)  to  (87)  and  (88)  and  the  radial  profiles  obtained  as  outlined 
in  part  A of  this  section. 


An  alternate  form  of  the  Abel  inversion  is 


(92) 


Although  apparently  simpler  than  the  form  of  (90)  and  (91),  the  latter  i» 
used  because  it  expresses  the  integration  portion  (91)  of  the  Abel  transform 
in  exactly  the  same  form  as  all  the  previous  integrals  that  have  been  con- 
sidered and  thus  the  same  quadrature  approximation  (see  part  E of  this 
section)  can  be  employed.  In  application  the  differentiation  portion  (90) 
of  the  Abel  transform  is  handled  with  a central  difference  approximation 
except  at  r = R where  a backcvard  difference  Is  uced,  and  at  r = 0 where  a 
symmetric  (about  r = 0)  quadratic  fit  to  F(r)  values  in  the  vicinity  of  r is 
used  to  obtain  g(r)  at  r = 0. 

The  Abel  inversion  xnetbod  can  be  used  to  o'ntain  meaningful  re  units 
only  for  optically  thin  sources.  Generally,  this  means  that  the  transverse 
txansmittance  profile  t(z)  be  greater  than  ''’0.90  for  ail  z.  For  smallsr 
ti ansmittance,  sortie  iterative  procedures  must  be  employed.  One  such 
procedure  La  considered  next  in  this  section.  The  Abel  inversion,  however. 
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plays  two  u'nportant  roles  in  this  iterative  method;  it  is  used  to  start  the 
iteration  process  and  to  continue  the  process  by  performing  inversions  of 
"difference”  profiles  obtained  by  subtracting  intermediate  calculated 
transverse  profiles  from  the  innut  transverse  profiles. 


C,  Iterative  Abel  Inversion 

In  order  to  obtain  inversions  for  radiation  sources  that  are  not 
optically  thin,  it  is  necessary  to  use  an  iteration  scheme.  The  method 
presented  here  was  inspired  by  a scheme  devised  by  Elder  et  al.  ' for 
the  inversion  of  monochromatic  emission/ absorption  data,  but  is  substan- 
tially modified  to  apply  to  band  n.odel  radiation  formulations  and  to  pro- 
vide a rational  test  for  convergence.  (That  the  formulation  used  by  Elder 
et  al.  can  give  convergence  to  the  correct  answer  is  highly  suspect.) 

The  starting  points  for  the  method  are  the  radiative  transfer  equa- 
tions (81)  and  (82).  Suppose  that  we  have  some  -'stimate  of  the  radial 
functions  c(r),  p{r),  and  T(r).  Denote  these  profiles  and  any  functions 
computed  from  them  by  the  subscript  i.  Then  we  can  write  from  (81)  and 
(82)  that 


R 


N (2)  = 2 \ J (r)wf(2.r)  - 

' I I ^“2  2 

z ''  r - z 


(93) 


R 


W.(z)  = 2 ^ K.(r)y;'(z,r)  — 

„ V 


rdr 


(94) 


2 2 
r - 2 


The  difference  profiles  obtained  by  subtracting  N.(z)  and  W^(2)  from  the 
known  profiles  N(2)  and  W(z)  are 

R 

an. (2)  = N(z)  - N.(z)  = 2 ^ 


J(r)w®(z,r)  - J.^(r)w?(z,  r) 


rdr 


/2  2 

''  r • z 


(95) 


iW.^(z)  = W(z)  - W.^(z) 


R 


z 


K(r)y®(z,r)  - K^(r)y^(z,r) 


rdr 


(96) 


wh^re  J(r),  K(r),  w®(z,r),  and  y^(z,  r)  are  functions  of  the  true  radial  pro- 
files c(r),  p(r),  and  T(r).  The  standard  procedure  for  obtaining  an  iteration 
formulation  is  to  assume  that  the  true  functions  will  obtain  in  the  next 
iteration.  Thus,  we  add  the  subscript  i + 1 to  all  nonsubscripted  variables 
on  the  right-hand  side  of  (95)  and  (96)  to  obtain 


AN.(z) 


R 
2 \ 


J.^j(r)w?^j(z,  r)  - J^(r)w®(z,r) 


rdr 


(97) 


AW.(z)  = 


R 

5 


rdr 


vTl 

j r 


(98) 


Despite  the  fact  that  the  coordinate  z appears  in  the  integrands  of  (97)  and 
(98),  we  now  apply  an  Abel  transformtition  to  ANj^(z)  and  iiW.^(z)  and  denote 
the  results  by  AJ|(r)  and  AK(r),  respectively.  In  some  sense,  we  then 
have 


or 


AJ. (r)  = J.^j(x->iwj^,  .(z,  r)  - J^(r)w®(z,  r) 

AK. (r)  = K.^j(r)y?^j(z,  r)  - K.(r)y?(z,r) 


J. 


i+i 


K 


i+1 


(r) 

(r) 


w.^j(z,r) 


yi+i(^,r) 


J. (r)w?(z,  r)  + 6J^(r/ 

K. {r)y°(z,i)  + &K.(r) 


(99) 

(100) 


(101) 


(102) 


Equations  (101)  and  (102)  define  an  iteration  procedure  for  J and  K,  but 
not  yet  in  a usable  form.  In  order  to  extract  a usable  form,  we  assume 
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a priori  tha^-  the  algorithm  converges.  Then,  as  the  iteration  progresses 

0_i.0  j0__..0  ,T,, 

w.,, — ►w.  andy.., — *"7..  Thus 
i+l  i •'i+l 


^i+l  = JL(r)  ^ AJ.(r) 

w.^j(z.r) 


(103) 


= K .:r,  *- 


(z»  r) 


AK.(r) 


(104) 


Further,  since  we  have  assumed  that  the  algorithm  will  converge,  we  have 

that  /iJ.(r)  — *-0  and  AK.(r) — ►0.  Thus,  it  will  make  no  difference  in  the 

^ ^ I'  0 

limit  whether  we  correct  the  i—  Iteration  for  J and  K by  AJ./w,  . , and 
0 i i+l 

^K./y.  . or  by  these  correction  values  multiplied  by  some  constant. 

^ 0 0 

Choosing  and  respectively,  as  the  multiplying  factors  for  the 

J and  K increments,  we  obtain  the  final  iteration  formula  as 

J.^j{r)  = J.(r)  + M.(r)  (105) 


Ki^l(r)  = K.(r)  + iK.(r) 


(106) 


where  AJ.(r)  is  the  Abel  transform  of  N(z)  - N.(z),  AK^(r)  is  the  Abel 
transform  of  W(z)  - Wj^(z);  Nj^(z)  is  obtained  from  (93),  and  W.(z)  is  obtained 
from  (94).  The  iteration  is  initiated  by  obtaining  "first  guess”  profiles  of 
c(r)  and  T(r)  by  Abel  inversion.  The  iteration  is  continued  until  some  con- 
vergence criteria  is  met.  Generally,  the  criterion  is  used  that  the  i—  and 
st 

i + 1 — radial  profiles  differ  by  no  more  than  soir.e  preset,  small  incre- 
ment. 

Although  this  iteration  scheme  has  been  derived  heuristically,  it  is 
shown  by  example  in  Section  VII  that  the  scheme  is  rapidly  convergent  and 
unique.  A summary  of  the  method  is  presented  as  a computational  flow 
diagram  in  Fig.  2.  Note  that  in  the  scheme,  the  pressure  profile  is  entered 
as  a known  profile. 

One  aspect  of  the  scheme  diagramed  in  the  figure  has  not  yet  been 
discussed  — the  testing  of  J(r)  and  K(r)  for  nonnegative  values.  If  tlie  input 


Fig.  2.  Iterative  Aoel  Inversion  Flow  Diagram 
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f 

r 


data  profiles  N (z)  and  t(z)  are  noisy  (i.  e.  , random  noise  fluctuations),  it 
is  possible  for  intermediate  results  for  the  iterated  J(r)  and  K(r)  profiles 
to  be  negative  If  they  are  so,  the  inversion  of  J(r)  and  K(r)  for  c(r)  and 
T(r)  is  susijcct  if  not  indeterminate.  In  the  inversion  scheme,  therefore, 
J(r)  and  K(r)  are  tested  at  each  discrete  value  of  r used  in  the  numerical 
approximation  routines  and,  if  found  to  be  nonpositive,  replaced  by  their 
value  from  the  preceding  iteration.  This  procedure  cannot  be  used,  of 
course,  if  negative  or  zero  values  of  J(r)  or  K(r)  occur  in  the  initiating 
Abel  inversion.  In  this  case,  the  input  transverse  pro^.l^s  must  be 
smoothed  to  such  an  extent  that  this  condition  does  not  arise.  Smoothing 
of  the  input  data  is  the  subject  of  the  following  part  of  this  section. 

D.  Pre  smoothing 

The  integral  equations  (81)  and  (82)  can  be  classed  as  Fredholm 
integral  equations  of  the  first  kind.  The  solution  of  such  equations  is 
known  to  be  highly  sensitive  to  irregularities  in  the  input  functions  N(z) 
and  W(z).  Even  very  small  fluctuations  (such  as  random  experimental 
noise)  in  the  input  functions  can  be  dramatically  amplified  in  the  inversion 
process  and  yield  essentially  meaningless  results.  Additionally,  there  is 
the  possibility  that  an  inversion  cannot  even  be  effected  if  the  Juctuaticns 
are  large  enough  to  cause  negative  or  zero  values  of  the  radial  functions 
J(r)  and  K(r),  Several  methods  have  been  devised  to  handle  this  fluctua- 
tion problem.  Here  we  consider  a straightforward  presmoothing  of  the 
input  data  to  be  inserted  into  the  Abel  or  iterative  Abel  inversion  routines. 

The  presmoothing  scheme  used  here  is  inspired  by  the  inversion 
schemes  of  Phillips^^^^  and  Twomey.  Consider  a function  g(z)  defined 

on  an  equal  interval  grid  z^  = (n  - l)Az  (n  = 1,  2,  . . . , N + 1)  consisting  of 
N intervals  of  width  Az.  Define  g^  = The  array  represents  either 

of  the  unsmoothed  transverse  profiles  K(z)  or  W(z),  and  Az  = R/N.  Let 
the  array  f^  (n  = 1,  2,  . . . , N + 1)  represent  the  smoothed  version  of  g^. 
The  Phillips -Twomey  approach  to  smoothing  is  to  select  that  set  of  f^ 
which  minimizes  the  curvature  function 
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N+1 


• • • ’ Wi^ 


E(f  . -Zi 

n-1  n 


+ f ) 
n+r 


(107) 


i=l 


subject  to  certain  constraints.  Clearly,  some  restraint  must  be  imposed 

since  the  set  that  minimizes  C without  constraint  is  the  trivial  set  f = 

n 

Constant  for  all  n.  The  constraint  imposed  here  is  that  the  root-mean- 
square  (rms)  difference  between  the  smoothed  and  unsmoothed  functions 
be  of  the  order  of  the  estimated  experimental  rms  error  e.  We  will  write 
this  constraint  as 


S(li.  • • • 


N+1 

(£^  - g/ - (N  t 1)^2 

i=l 


(108) 


The  method  of  Lagrange  undetermined  multiplier  indicates  that  we  should 
then  minimize  the  function  C - XS  with  respect  to  the  f^  and  then  vary  X 
until  the  constraint  (108)  is  satisfied.  Rather  than  X,  we  use  v = -1/X  and 
minimize 


1 ^N+1^  "" 

^^^^1 ^N+l^  • • »^N+l' 

N+1 

N+1 

= ''  2 ‘'n-l  - 

2 2 

n n+ 1 n n 

(109) 

n=  1 

n=l 

Minimization  of  H with  respect  to  the  f^  is  obtained  by  solving  the  system 
of  N + 1 equations 


bH 

bf, 


= 0 


bH 

bf. 


= 0 


(110) 


bH 


bf, 


N+1 
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Differentiation  of  (109)  with  respect  to  L,  use  of 

bf 

" - A 

bf.  nj 

J 


(111) 


and  equating  the  result  to  zero  gives 

^ ^ ^ ^V2  ~ 

This  result  is  the  form  of  the  equations  (110)  .or  j = 3,  4,  . , N - 2. 

For  j in  the  end  regions  j = 1,  2 and  j = N-  l,  N,  N+1,  the  form  of  the 
equations  will  depend  on  how  f^  is  defined  for  n ^ 0 and  n ^ N + 2.  The 
transverse  profiles  N(z)  and  W(z)  are  symmetric  about  z = 0 and  fall 
(ideally)  to  zero  at  z = R.  It  seems  appropriate  to  define 

f.  = f2_.  j . 0 (113) 


^N+1  " 


0 


= ■^2(N+l)-j 


j a N + 2 


(114) 


A common  approach  is  to  substitute  the  conventions  directly  into  (112)  to 
obtain  the  special  end  zone  equations.  However,  this  procedure  is  not 
strictly  correct  because  then  we  are  assuming  a dependence  between  f^ 
values  after  we  have  already  used  the  independence  relation  (111).  The 
correct  procedure  is  to  use  (113)  and  (114)  in  H before  the  partial  differ- 
entiations are  taken.  Thus, 


H = 


(-2£,  + 


N-1 


(f 


n-1 


- 2f  + 


n 


n=2 


N 


('n- 


'n 


.2  ^ 2 
^ ®N+1 


(115) 


From  this  expression  for  H,  the  following  system  equations  result 
j = 1 (5v  + l)fj  - = gj 

j = 2 -6Yfj  + (9y  + 1)^2  ■ 

(116) 

J " 1 Yf^.s  - 4Yfj^-.2  + (6y  + - 4Yf^  = 

j = N Vfj^_2  - 4Yfj^_j  + (5y  + l)fj^  = gj^ 

Note  that  because  been  specifically  set  to  zero,  the  system  con- 

tains only  N equations.  With  (112)  and  (-16),  the  system  of  equations  can 
be  written  in  the  matrix  form 

(yH  + T)?  = i (117) 

where  f and  g are  the  column  matrices 


and  H ir,  the  banded,  N X N,  symmetric  matrix 


5 

-6 

1 


H = 


-6 

9 

-4 

1 


1 

-4 

6 

-4 


1 

-4  1 

6 -4 


1 


1-4  6-4  1 

1-4  6-4 

1 -4  5 


The  formal  solution  for  the  smoothed  array  f^  is  then,  from  (117) 


f = (yH  + I)-^  g 


(120) 


(121) 


The  smoothing  parameter  y controls  the  degree  of  smoothing.  For  y = 0, 

f = I ^g  = Ig  = g and  no  smoothing  is  effected.  As  y — >■  ®,  f — ►y  ^ g- — ►O. 

In  practice,  y is  varied  until  the  constraint  (108)  is  met.  Application  of 

this  smoothing  method  requires  a computation  of  the  inverse  of  yH  + I. 

(22) 

A standard  inversion  routine  based  on  the  Gauss  elimination  algoritnm 
is  employed.  Characteristics  of  this  smoothing  routine  in  application  to 
real  data  are  discussed  in  Section  VIIIE. 


For  some  applications,  it  is  necessary  to  smooth  data  defined  on  an 
irregular  grid.  The  generalization  of  the  method  considered  above  to  a 
general  grid  is  straightforward,  but  tedious.  Only  an  outline  of  the  deriva- 
tion and  the  final  result  are  reported  here.  The  unsmoothed  (and  smoothed) 
profiles  are  taken  to  be  defined  on  the  general  grid  z^,  i = 1,  N + 1 with 
z = 0 and  - R-  Define  the  interval  weighting  functions 

^i  ^ *i+l  ■ ^i-1  i=2,  3,  ...,N  (122) 

‘^N+l  ^ ^N+1  ■ 
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The  generalization  of  (109)  for  the  function  to  be  minimized  is 


N+1 


H(f 


f»,, ,)  = Y f 1 " b f + c f . i]^ 

1 N+1  / j - n n-1  n n n n+l-* 


n=l 

N+1 


(f  - g )^d 
/ j n “n  n 


n=l 


(123) 


where  the  second-derivative  coefficients  a , b and  c are 

n n n 


n ~ (z  , . - z ,)(z  - z ,) 

' n+1  n-1  n n-r 


^n  ~ (z  , , - z )(z  - z ,) 

n+1  n n n-1' 


y n = 1,  2,  . . , , N 


(124) 


with  the  convention  that  z^  s -z^  and  = 2R  - Zj^j.  Carrying  the  deriva- 

tion through  as  before,  we  find  the  smoothed  profile  to  be  given  by 


f = K'^  g' 


(125) 


where  g’  is  the  column  matrix 


g'  = 


and  K is  the  banded,  N x N,  symmetric  matrix 


‘iig 


1^1 


*^2^2 


(126) 
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with 


E.  Quadrature  ApproxlTnatlon 

All  of  the  integrals  required  in  radiation  computation  and  inversion 
can  be  written  in  the  general  form 


F(y,  P) 


*\  G(y,  x) 


2 2 
X - y 


(129) 


y P s:  R 


(130) 


For  integrals  over  only  a portion  of  an  optical  path  [such  as  (19)  through 
(21)  and  (28)  to  compute  path  averaged  band  model  parameters],  we  set 
X = r,  y = z and  P = r.  For  integrals  along  an  entire  optical  path  [such 
as  (6)  and  (7)  v;hich  compute  the  radiance  and  transmittance  from  the 
chord  at  transverse  coordinate  z],  we  set  x = r,  y = z and  P = z.  For  the 
Abel  inversion,  we  set  x = z,  y = r and  P = r. 


The  quadrature  approximation  for  (129)  is  patterned  on  the  method 

(23) 

■r.  Divide  the  interval  x = 0 to  x = R into  N equal  size  subinter- 


of Barr. 


vals.  The  subinterval  size  is  then 


A R 
^ ~ N 

The  coordinate  x and  parameters  y and  P are  discretized  as 


x^  = (n  - 1)A 


y,  = (k  - 1)A 


n = 1,  2,  . . . , N + 1 


k=l,  2,  N+1 


(131) 


(132) 


Pq  = (q  - DA 


q=l,  2,  N+1 


where,  because  of  (130)  and  the  integration  lunits  on  (129) 

Isk^qSnsNr 1 


(133) 


Define  the  discrete  function  values 


F,  = F(y,  , P ) 
kq  '^k  q' 


(134) 


Within  the  subintervai  ^ x S the  kernad  function  G(y,x)  is  approx- 

imated by  the  symmetric  quadratic  form 


G(y.x)  = + 


(135) 


The  coefficients  aj_^  and  are  determined  by  forcing  G(y,x)  to  be  equal 

to  G,  and  G,  , , at  the  interval  end  points  x and  x . respectively, 
kn  k,  n+  i ^ n n+ 1 ^ ^ 


The  solutions  are 


kn 


^ (2n  - 1) 


'«kn  - ” ^)'^k 


,n+lj 


A^(2n  - 1) 


n+ 1 " ^kn 


(136) 


With  the  relationships  so  far  established,  (129)  is  written  as 


F,  =2 

kq 


N nC^ 
r 


Z \ 

n=q  (n-l)A 


a,  + b,  X 
Kn  kn 


xdx 


/~2  ? 
^ - ^k 


(137) 


Evaluation  of  the  integrals  over  x and  a straightforward  but  tedious  manip- 
ulation of  the  resulting  expressions  gives  the  final  result 

N 

Fr,  = (1  <>j(k,q)GRq+  2 ''2<‘=-''l°kn  + ‘'3*^>°k.N+l 

L n=qT 1 

where  the  quadrature  coefficients  are 


(133) 
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ajCk.q)  = 3 


2 2 
q - (k-ir 

2q  - 1 


■IV'"' 


"ri 


1)^  - (k-i)^  - 


2^/2 


2a  - 1 


(139) 


2 2 I 2 2 

n - (k-1)  _ 4(n-l)  l(n-l)^  - (k-1)^ 


13/2 


2n  - 1 


(2n-l)(2n-:^) 


(n-2)^  - (k-1)^ 


3/2 


2n  - 3 


(140) 


“ 3 


!(N-1)^  - , 3 / 2 2*  iN^-(k-l)^| 

2N  - 1 2 V ^ ' 2N  - 1 


3/2 


(141) 


The  only  qualifications  required  for  application  of  (138)  are  the  index 
relation  of  (133),  the  convention  that  the  summation  term  be  omitted  if 
q = N and  the  obvious  result  that  = 0ifkorq  = N+  1. 
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IV.  PROGRAM  EMABIC 


A.  Description  of  Code 

A computer  code,  EMABIC,  has  been  written  which  incorporates 
the  band  model  radiation  formulation  of  Section  II  and  the  inversion  and 
smoothing  formulations  of  Section  III  into  a general  usage  program.  The 
code  is  written  in  Fortran  IV  compatible  with  The  Aerospace  Corporation 
CDC  7600  and  the  AFRPL  CDC  6400  computer  systems. 

The  code  operates  in  one  oi  four  modes,  a transverse  profile  gen- 
eration mode,  a data  smoothing  mode,  and  two  data  inversion  modes. 

The  specification  of  mode  and  data  to  be  analyzed  is  treated  in  part  B of 
this  section.  Here,  a brief  description  of  the  operation  of  the  code  in 
each  of  these  modes  is  given. 

A routine  common  to  all  modes  is  ti...  tput  of  relevant  data.  This 
input  is  accomplished  by  calls  to  the  two  subroutines  INPUT  and  DAT- 
PREP.  The  first  subroutine  simply  reads  in  all  data  required  for  a com- 
putational run,  stores  the  data,  and  provides  a preliminary  output  listing 
summarizing  the  calculation  conditions  and  data.  The  second  subroutine 
prepares  the  input  data  for  use  by  the  main  program  EMABIC,  The  two 
most  important  preparations  are  the  fitting  of  input  protiles  to  an  equal- 
interval  grid  on  the  coordinate  r or  z in  the  range  zero  to  R and  smooth- 
ing of  input  transverse  data  profiles.  The  former  is  accomplished  by  a 
call  to  subroutine  ZONEFIT  which  fits  the  input  profile  (defined  on  an 
arbitrary  coordinate  grid)  by  linear  interpolation  to  the  grid  (i  - 1)  R/N, 
i = 1,  2,  . . .,  N + 1 where  N is  the  number  of  zones  used  in  the  computa- 
tions. Smoothing  is  performed  by  a call  to  subroutine  SMOOTH  which 


EMABIC  is  an  acronym  for  Emission/Absorption  Inversion  Code. 

Requests  for  this  code  by  interested  users  should  be  directed  to  AFRP.L/CE 
in  writing. 
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employs  the  formalism  of  Section  IIID.  The  matrix  inversion  of  (121)  is 
handled  by  subroutine  GELB.  Provisions  are  also  made  in  subroutine 
DAT  PREP  for  the  removal  of  finite  field-of-view  effect'i  on  transverse 
profiles  by  a call  to  subroutine  DCONVLV, 

A second  conimon  routine  is  the  computation  of  the  quadrature 
weighting  coefficients  (139-140)  by  a call  to  subroutine  COEFF,  These 
coefficients  are  subsequently  used  in  any  call  to  subroutine  QUAD  which 
applies  the  quadrature  approximation  of  (138)  whenever  a numerical 
integration  is  called  for. 

Profile  Generation  Mode.  The  profile  generation  mode  accepts  input 
radial  profiles  of  pressure,  temperature  and  concentration  and  computes 
transverse  profiles  of  radiance  and  transmittance  according  to  the  band 
model  radiation  formulation  of  Section  U.  Both  N{2)  and  t(2)  are 
computed  by  a single  call  to  subroutine  PROFILS.  Depending  on  the  approxi- 
mation used  to  treat  optical  path  nonuniformities  and  the  spectral  lineshape, 
one  or  more  of  the  following  routines  may  be  called  by  PROFILS. 

YCGL  Compute  the  derivative  function  y for  the  CG  approximation 
and  Lorentz  lineshape. 

YDRL  Compute  the  derivative  function  y for  the  DR  approximation 
and  Lorentz  lineshape.  YDRL  calls  the  subprogram  YLSL 
which  computes  y for  the  LS  approximation  and  Lorentz 
lineshape. 


^Subroutines  GELB  and  GELS  are  taken  from  the  IBM  Scientific  Sub- 
routine Package  (Reference  22). 

DCONVLV  has  not  yet  been  coded  and  at  present  an  immediate  return 
to  the  main  program  takes  place  in  a call  to  DCONVLV. 
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YCGD  Compute  the  derivative  function  y for  the  CG 
approximation  and  Doppler  lineshape. 

YMLD  Compute  the  derivative  function  y for  the  ML  derivative 
approximation  and  Doppler  lineshape.  YMLD  calls  the 
subprogram  YLSD  which  computes  y for  the  LS  approxi- 
mation and  Doppler  lineshape. 

YMDC  Compute  the  derivative  function  y for  the  Voigt  lineshape. 

F Compute  the  curve- of-growth  function  for  a band  of 

Lorentz  lines 

G Compute  the  curve-of-growth  function  for  a band  of 

Doppler  lines 

WMIX  Compute  the  curve-of-growth  function  for  a band  of 
Voigt  lines. 

Band  model  parameters  required  for  generation  of  transverse  profiles 
are  obtained  from  a call  to  subroutine  KDPARAM  which  interpolates 
linearly  in  temperature  on  a table  of  parameters  entered  as  part  of  the 
input.  The  results  of  a profile  generation  run  are  listed  by  subroutine 
LIST  I . 

Profile  Smoothing  Mode.  For  an  inversion  run  using  the  Abel  or  iterative 
Abel  method,  the  input  transverse  profile  data  must  first  be  smoothed.  For 
the  actual  inversion  run,  this  smoothing  is  part  of  the  data  preparation 
performed  by  subroutine  DATPREP  and  uses  smoothing  parameters 
entered  as  part  of  the  input.  In  order  to  determine  the  proper  smoothing 
parameters  for  the  inversion,  a prior  run  of  EMABIC  in  the  variable  smooth- 
ing mode  is  generally  required.  In  this  mode,  the  transverse  profiles  are 
s moothed  according  to  the  procedure  of  Section  UID  for  a wide  range  of 
smoothing  parameter  "Y  . Inspection  of  the  degree  of  smoothing  performed  as 
a function  of  V allows  the  selection  of  proper  V values  for  the  actual 
inversion.  The  principle  subroutine  used  in  this  mode  is  subroutine  GAMMA 
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Abel  Inversion  Mode.  When  operating  in  this  mode,  EMABIC  inverts  the 
input  transverse  profiles  N(2)  and  t (z)  for  c(r)  and  T(r)  according  to  the 
Abel  inversion  routine  of  Section  II^B.  The  first  step  of  the  inversion  is 
the  computation  of  the  radial  source  and  absorption  functions  J(r)  and  K(r). 
These  calculations  are  effected  by  a call  to  subroutine  ABEL.  The  radial 
profiles  T(r)  and  c(r)  are  obtained  from  J(r)  and  K(r)  according  to  (85)  and 
(86)  in  subroutine  CTSOLVE.  K the  solution  for  c{r)  and  T{r)  is  indeter- 
minate because  J(r)  is  negative,  an  error  message  is  written.  Generally, 
this  indeterminacy  can  be  removed  by  a higher  degree  of  data  smoothing. 

The  results  of  an  Abel  inversion  are  written  out  by  subroutine  LIST3. 

Iterative  Abel  Inversion  Mode.  In  this  mode,  inversion  of  transverse 
profiles  is  carried  out  with  the  iterative  Abel  inversion  method  derived  in 
Section  lUC.  The  successive  Abel  inversions  on  the  difference  profiles 
and  profile  generations  using  updated  radial  profiles  (see  Fig.  2 ) are 
carried  out  with  subroutines  ABEL  and  PROFir,S,  respectively.  At  each 
iteration,  the  source  functions  J(r)  and  K(r)  are  used  to  obtain  c(r)  and 
T(r)  as  in  the  Abel  inversion.  Testing  for  convergence  is  made  in  subroutine 
CONVERGE.  Convergence  criteria  are  set  as  part  of  the  input  data  and 
provide  for  testing  on  the  maximum  absolute  difference  or  root-mean- 
square  difference  between  successive  iterations  of  any  or  all  of  the  four 
profiles  c(r),  T(r),  N(z)  or  t(z).  Results  of  a successful  inversion  (and, 
optionally,  the  results  of  intermediate  iterations)  are  written  by  subroutine 
LIST2.  Optionally,  intermediate  iteration  results  for  J(r)  and  K(r)  may  be 
listed  by  subroutine  LIST4. 
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B.  Preparation  of  Input  Data 

A computational  run  of  program  EMABIC  requires  a s et  of  program 
control  cards  to  specify  the  mode  of  computation  and  to  supply  input  data. 
Some  program  control  cards  simply  specify  a computation  mode,  some 
specify  a computation  mode  and  supply  data,  while  others  signal  the  code 
that  blocks  of  auxiliary  data  are  now  to  be  read  in.  Each  type  of  control 
card  contains  an  alphanumeric  name  in  the  first  ten  card  columns.  These 
names  must  be  spelled  correctly  and  must  be  left- justified.  If  data  are 
specified  on  a program  control  card,  they  must  be  entered  in  accordance 
with  the  format  specification  indicated  in  the  detailed  description  of  each 
card  given  below.  All  fields  of  the  program  control  cards  are  10  columns 
wide.  In  general,  integer  and  alphanumeric  data  must  be  right-justified 
in  their  fields.  Non-integer  numerical  data  may  be  entered  in  either  F or 
E formats  (with  decimal  point  and,  for  the  latter,  the  exponential  symbol  E). 
E-formated  data  must  be  right- justified  in  their  field.  These  same  rules 
apply  to  data  entered  on  auxiliary  card  decks.  The  types  of  control  cards 
and  the  data  contained  on  them  are  illustrated  in  Fig.  3.  A description  of 
each  type  follows. 

1.  Title  Card.  The  card  name  is  TITLE.  Columns  11-80  of  this  card 
may  be  used  for  any  identification  title  desired.  This  title  will  appear  on 
the  final  listing  of  the  computation  run. 

Z.  Lineshape  Card.  The  card  name  is  LINESHAPE.  The  variable 
SHAPE  (format  AlO)  defines  the  lineshape  that  will  be  us  ed  in  the  statistical 
band  model  radiation  calculations  and  must  have  one  of  the  three  alpha- 
numeric values  LOP.ENTZ,  DOPPLER  or  VOIGT. 

3.  Nonuniformity  Approximation  Card.  The  card  name  is  INHOM  (for 
inhomogeneous).  The  variable  APPROX  (format  AlO)  specifies  the 
approximation  that  will  be  used  to  treat  nonuniformities  along  an  optical 
line  of  sight.  If  APPROX  has  the  value  CG,  the  Curtis-Godson  approxima- 
tijn  will  be  used.  If  the  value  is  DR,  the  derivative  approximation  will  be 
used. 
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4.  Integration  Grid  Card.  The  card  name  is  ZONES.  The  variable 
NZONES  (format  110)  defines  the  number  of  equal  size  radial  and  transverse 
zones  that  will  be  used  in  the  various  numerical  routines  of  EMABIC.  The 
maximum  allowed  value  is  50. 

5.  Band  Model  Parameter  Card.  The  card  name  is  PARAM.  This  card 
calls  for  the  read-in  of  band  model  parameters  for  the  active  species  and 
spectral  bandpass  being  considered.  The  auxiliary  data  defining  the  band 
model  parameters  must  follow  the  PARAM  card  immediately  in  the  card 
input  stream  and  must  be  structured  as  shown  in  Fig.  4.  If  the  variable 
PRINT  (format  AlO)  has  the  value  PRINT,  the  band  model  parameters  will 
be  listed. 

6.  Transverse  Profile  Mode  Card.  The  card  name  is  PROFILES. 

This  card  is  used  to  operate  EMABIC  in  the  profile  generation  mode.  That 
is  , transverse  profiles  of  radiance  and  transmittance  will  be  generated 
from  input  radial  profiles  of  pressure,  temperature  and  species  concentra- 
tion. This  input  data  mus  t immediately  follow  the  PARAM  card  and  be 
structured  as  illustrated  in  Fig.  5.  If  the  variable  PRINT  (format  AlO)  on 
the  PROFILES  card  has  the  value  PRINT,  these  radial  data  will  be  listed. 
The  auxiliary  radial  data  deck  allows  for  the  specificatior  of  profiles  for 
H2O,  CO^  and  CO.  In  order  for  the  program  to  know  which  species  to 
work  with,  a SPECIE  control  card  (type  7)  must  be  used  with  a PROFILES 
card. 

7.  Species  Card.  The  card  name  is  SPECIE.  This  card  must  be 
included  as  part  of  the  control  card  sequence  whenever  EMABIC  is 
operating  in  the  transverse  profile  generation  mode.  The  variable  SPECIE 
(format  AlO)  must  have  one  of  the  three  alphanumeric  values  H20,  COE  or 
CO.  Note  that  this  control  card  only  tells  the  program  where  to  find  the 
correct  concentration  data  in  the  radial  input  data  deck.  If  does  not  set 
band  model  parameters,  for  example.  Thus,  the  band  model  parameters 
entered  with  a PARAM  card  (type  5)  must  be  consistent  with  the  species 
called  for  by  the  SPECIE  card. 
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Legend 


All  field  formats  are  ElO  or  FIO  except  the  PARAMID  field 
which  IS  AlO. 


PARAMID 

V 

Av 

■7 


“1 


"2 

"3 

T(i) 

k{  i ) 

D(  i ) 


Any  lO-character  identification  name. 

Spectral  position  (cm  ). 

Spectral  resolution  (cm  S-- 

Pressure  broadening  coefficient  (cm  * 'atm  for 
non-resonant  self-broadening  at  STP)., 

Ratio  of  resonant  self- broadening  parameter  10  Y 
at  STP. 

Ratio  of  foreign-gas  broadening  parameter  to  Y 
at  STP. 

Atomic  weight  of  active  gas  species  (amu  ). 
Temperature  array  ('^K).  The  T(  i ) array  must 

be  T(  i ) = lOOi,  1 = 1,2 40. 

Absorption  coefficient  for  v , A."  and  T(  i ) 

(cm  */atm). 

Line  density  parameter  tor  v , Av  and  '!(  1 ) 
(lines/cm  S..  Note,  D = 1/6  . 
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Fig.  4.  Input  Card  File  Structure  for  Band  Model  Parameters 


-60- 


legend 


All  field  formats  att:  ElC  or  FIO  except  the  DATAiU  field  which 
AlO  and  the  NPOlNTS  field  which  is  liO. 


lU 


DATAID 

NPOlNTS 


R 

r(  1 } 


pi  1 ) 
Tl  1 } 
^HZO 
^COZ 

cco 

^0 


Any  lU-character  identification  name. 

N.  mber  of  radial  positions  at  which  data  are 

given  (ma.*c2mum  value  = lOl), 

Source  cylinder  radius  (cm). 

Radiol  positions  at  which  data  are  given  (cm). 
The  value  of  r(  1)  must  be  rero,  the  value  of 
r(NPOiNTS)  must  be  R,  and  r(  j ) must  increase 
w ith  me  rcasipg  i . 

Pressure  at  r(  i l (atm). 

Temperature  at  r(  r ) (*^K). 

Concentration  of  H2O  at  r(  1 ) (mole  traction) 
Concentration  of  CO^  at  r(  . ) (mole  fraction). 
Concentration  of  CO  at  r(  1 ) (mole  fraction). 
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Fig.  5.  Input  Card  File  Structure  for  Radial  Profile  Data 


8.  Inversion  Mode  Card.  The  card  nanne  is  INVERSION.  An  encounter 
of  this  card  in  the  input  stream  signals  EMABIC  that  an  inversion  run  is 
to  be  made.  Usually,  a card  deck  of  data  defining  the  transverse  radiance 
and  transmittance  profiles  that  are  to  be  inverted  must  follow  this  card 
immediately  in  the  input  stream.  The  deck  structure  for  these  data  is 
shown  in  Fig.  6.  If  the  variable  PRINT  'format  AlO)  on  the  INVERSION 
card  has  the  value  PRINT,  these  input  data  will  be  listed.  The  exception 
to  the  requirement  that  an  auxiliary  data  deck  immediately  follow  the 
INVERSION  card  occurs  if  the  variable  SAME  (format  AlO)  has  the  value 
SAME.  In  this  case,  the  transverse  profiles  data  that  were  used  in  the 
immediately  preceding  run  will  be  used  unchanged.  The  variable  MODE 
(format  AlO)  specifies  which  of  two  inversion  methods  will  be  used.  If 
MODE  has  the  value  ABEL,  the  optically  thin  Abel  inversion  will  be 
employed.  If  the  value  is  INVERSION,  the  iterative  Abel  method  will  be 
used.  The  variables  G1  and  G2  (each  format  ElO)  are  smoothing  para- 
meters for  the  absorption  and  radiance  input  profiles,  respectively.  For 
an  iterative  Abel  inversion,  these  smoothing  parameters  are  presmoothing 
parameters  and  are  determined  from  a prior  run  of  EMABIC  in  the 
smoothing  mode  (following  paragraph). 

A fourth  value  is  possible  for  the  variab3e  MODE.  If  the  value  is 
SMOOTH,  no  inversion  of  the  input  profiles  is  made,  but  rather  EMABIC 
performs  a variable  smoothing  of  the  input  transverse  radiance  and 
transmission  curves.  Both  profiles  are  smoothed  using  values  of  Y 
starting  at  G1  and  increasing  by  factors  of  two  up  to  and  including  G2. 

The  output  of  EMABIC  operating  in  this  mode  allows  the  selection  of  G1 
and  G2  for  operation  in  the  iterative  Abel  inversion  mode. 

9.  Pressure  Profile  Card.  The  card  name  is  PRESSURE.  An 
encounter  of  this  card  indicates  that  a card  deck  specifying  the  radial 
pres  sure  profile  that  will  be  used  in  an  inversion  computation  follows 
immediately  in  the  data  input  stream.  This  card,  and  the  data  it  calls 
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Legend 


All  field  formats  are  ElO  or  FIO  except  the  DATAID  field  which 
IS  AIO>,  the  NPOINTS  field  which  is  110,,  and  the  CZ  field  which 

IS  no. 


DATAID 

NPOINTS 


R 

A 2 


Cl 


C2 


2(i) 


N(i) 


•(i) 


An/  lO^character  identification  name. 

Number  of  transverse  positions  at  which  data  are 

given  (maximum  value  = 101). 

Source  cylinder  radius  icm). 

Transverse  spatiil  resolution  of  data  (i.e. 
instrument  FOV  resolution)  (cm).*" 

A constant  that  converts  by  multiplication  the 

— 2-1 
N(  i ) emission  data  to  units  of  W/tm  -sr-cm 

(If  Cl  IS  tero  or  blank, > Cl  ~ x is  assumed). 

A flag  that  indicates  whether  the  ’’absorption" 

da'a  18  absorptance  or  transmittance.  Blank  or 

zero  (normal  mode)  implies  transmittance. 

C2  - 1 implies  absorptance. 

Transverse  positions  it  which  data  are  given  (cn.). 

The  value  of  z(l)  must  he  zero,,  the  value  of  z(NPOlNTS) 
must  be  R + Az/^  and  z(  i 1 must  increase  with 
increasing  \ , 

Emission  data  at  zu).  (Nominal  unit  is  W/cm^*sr- 
\ see  use  of  Cl).  The  value  of  N(NPOIN’TS) 
must  be  zero. 

Absorption  data  at  z(i).  (Nominal  variable  is  trans- 
mitta-^ce,  see  us-*  of  C2),  The  value  Oi~{NPOF.TS) 
must  De  unity  if  the  data  are  transmittance  and  zero 
if  absorptance.  "fii)  must  not  be  either  zero  or 
uiiity  for  any  other  va’ue  of  i. 


Az  = 0 IS  the  inly  value  that  can  presently  be  used, 
1C  20  30  40 


50 


60 


SO 


DATAID 

NPOINTS 

R 

AZ  Cl  C2 

d(l) 

Nil) 

^(1) 



4(2) 

N(Z1 

: i 

N(NPOINTS) 

t(NPOINTS) 

Fig.  6.  Input  Card  File  Structure  for  Transverse  Profits  Data 
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for  must  be  included  whenever  an  inversion  mode  is  specified.  That  is,  when- 
ever an  INVERSION  card  ia  used  and  MODE  + SMOOTH.  The  card  file 
structure  for  the  pressure  data  is  given  in  Fig,  7.  If  the  variable  PRINT 
(format  AlO)  has  the  value  PRINT  on  the  PR^^SSURE  Control  card,  the 
radial  pressure  profile  will  be  listed. 

10.  Lis  t Control  Cards  . Two  cards  with  the  names  LIST  and  NO  LIST 

control  the  printout  of  transverse  and  radial  profiles  generated  during 
successive  iterations  in  an  inversion  mode  computation  If  a LIST 

card  is  never  encountered,  the  printout  of  intermediate  profiles  is  sup- 
pressed. If  a LIST  card  is  encountered,  the  intermediate  profiles  will  be 
listed.  A NOLIST  card  can  be  used  to  again  suppress  the  listing  on  a run 
subsequent  to  a run  in  which  the  listing  was  called  for.  If  the  variable 
JKLIST  (format  AlO)  on  the  LIST  card  has  tlie  value  LIST,  the  emission 
and  absorption  functions  J(r)  and  K(r)  are  also  listed  for  each  iteration. 

11,  Convergence  Card.  The  card  name  is  CONVERGE  and  contains 
data  specifying  ihe  convergence  criteria  that  must  be  met  in  order 

to  terminate  an  inversion  computation.  The  variable  CTEST  (format  AlO) 
must  have  one  of  the  two  values  TC  or  RT.  The  former  value  indicates 
that  convergence  testing  will  be  done  on  the  radial  temperature-concen- 
tration profiles  while  the  latter  indicates  that  the  tests  will  be  nr:ade  on 
the  transverse  radiance-transmittance  profiles. 

The  variable  IMAX  (format  XIO)  specifies  the  maximum  nvinber  of 
iterations  that  will  be  allowed.  If  convergence  is  not  achieved  within 
this  number  of  iterations,  a message  to  this  effect  is  printed,  the  radial 
and  tran.sverse  profiles  of  the  last  iteration  are  printed,  and  the  program 
moves  on  to  the  next  run. 

If  CTEST  has  the  value  TC,  the  variables  TRRMS,  CTRMS,  TRMAX 
and  CTMAX  (all  format  ElO)  have  the  following  meaning: 
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Legend 


All  fields  are  ElO  or  FIO  except  the  PPROFID  field  which  is 
AlO  and  the  NPOINTS  field  which  is  110. 


PPROFID 

NPOINTS 


r(  i ) 


p(  1 ) 


Any  10-character  identification  name. 

Number  of  radial  points  at  which  pres  sure  is 

defined  (maximum  value  = 101). 

Source  cylinder  radius  (cm). 

Radial  positions  at  which  pressures  are  given  (cm). 
The  value  of  r(l)  must  be  zero,  the  value  of 
r(NPOINTS)  must  be  R,  and  r(  i ) must  increase 
with  increasing  i . The  r(  i ) array  need  not 
(and,  in  fact,  cannot  unless  Az  = 0)  correspond 
with  the  transverse  input  data  grid  zC  i ). 

Pressure  at  r(  i ) (atm). 


Fig.  7.  Input  Card  File  Structure  for  Radial  Pressure  Profile 


The  absolute  value  of  the  maximum  difference  (at  any  radial  position) 
between  bvo  successive  temperature  profiles  must  be  less  than 
TRMAX  (unit  = K)  for  convergence  to  be  assumed.  The  absolute 
value  of  the  radially  computed  root  mean  square  difference  (rms) 
between  the  two  profiles  must  be  less  than  TRRMS  (unit  = K). 

Similarly,  the  maximum  absolute  difference  between  successive 
concentration  profiles  must  be  less  than  CTMAX  and  the  rms  dif- 
ference must  be  less  than  CTAVE. 

If  CTEST  has  the  value  RT,  TRRMS  ..r.J  TRMAX  correspond  to  the 
rms  and  maximum  absolute  differences,  r>.-i  nectively,  allowed  between 
two  successively  computed  transverse  radiance  profiles.  In  this  case, 
the  unit  of  TRRMS  and  TRMAX  is  W/cm^-sr-cm"  ^ CTRMS  and  CTMAX 
now  correspond  to  tests  mads  on  successively  computer  transmittance 
profiles. 

The  most  stringent  convergence  testing  occurs  when  all  four 
variables  (TRRMS,  TRMAX,  CTRMS,  and  CTMAX)  have  small  but  non- 
zero values.  A zero  value  for  any  of  the  variables  indicates  that  testing 
on  that  value  is  not  made.  If  all  four  variables  are  specified  to  be  zero, 
no  convergence  testing  is  made  and  the  results  are  assumed  to  converge 
after  the  first  iteration. 

12.  Execution  Card.  The  card  name  is  RUN.  When  this  card  is 
encountered,  computations  are  begun  using  the  data  entered  up  to  that 
point,  and  an  output  listing  of  the  results  is  made.  When  the  computation 
and  results  listing  is  completed,  the  program  continues  to  read  program 
control  cards  until  a new  RUN  card  is  encountered.  A new  computation 
is  then  begun  for  all  of  the  conditions  and  data  of  the  first  run  except  those 
which  have  been  changed  by  the  intervening  program  control  cards  and 
auxiliary  data  decks.  This  process  is  repeated  until  an  end-of-file  card 
is  encountered.  With  this  feature,  a large  number  of  related  runs  can  be 
made  with  one  job  submission. 


Other  than  the  requirement  that  all  required  data  be  specified  before 
a RUN  card  is  encountered  and  that  auxiliary  data  immediately  follow  the 
control  card  that  calls  for  them,  the  program  control  cards  may  be  arranged 
in  any  order. 

Great  care  should  be  taken  in  the  preparation  of  input  data  since 
very  few  checks  of  data  consistency  and  setting  of  default  values  are 
provided.  A general  feature  of  data  preparation  is  that,  if  particular 
data  on  a card  are  not  required,  they  need  not  be  specified.  If  none  of  the 
data  on  a control  card  is  needed,  that  card  need  not  be  included. 


! 


I 
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V.  BAND  MODEL  PARAMETERS 

Application  of  EMABIC  for  either  transverse  emission /absorption 
profile  generation  or  inversion  requires  the  specification  of  band  model 
parameters  for  the  active  species  of  interest.  The  basic  band  model  para- 
meters used  here  for  H^O  and  CO^  are  those  constructed  by  Young^^ 
explicitly  for  use  in  combustion  gas  problems  in  which  a high  degree  of 

nonisothermality  obtains.  Briefly,  the  parameters  were  constructed  by 

(181 

combining  the  high-temperature  parameters  from  the  NASA  Handbook  ' ’ 

with  low-temperature  parameters  derived  from  the  optical  line  data  of  the 
AFGL  compilation.'  ' The  former  parameters  are  quite  adequate  for 
temperatures  above  ~1500  K but  fail  for  some  lower  temperatures.  The 
parameters  derived  from  the  AFGL  compilation,  on  the  other  hand,  are 
very  good  below  ^500  K,  but  fail  at  some  higher  temperatures  because  of 
the  lack  of  hot-line  data  in  the  compilation.  The  combined  parameters 
represent  the  state-of-the-art  parameters  for  H2O  and  CO2  for  application 
to  the  temperature  range  200  to  3000  K. 

The  combined  parameters  k(T,v),  6(T,  v)  for  H2O  in  the  2.  T-pm  region 
are  tabulated  in  Ref.  23  and  24  for  T = 100,  200,  300,  500,  750,  1000,  1500, 
2000,  2500,  3000  K and  v = 2500  to  4500  cm"  ^ in  steps  of  25  cm~^.  The 
spectral  resolution  of  the  data  is  25  cm  The  same  parameters  for  CO2 
in  the  4.  3 -pm  region  are  tabulated  in  Ref.  25  for  the  same  temperature  array 
and  for  v = 2C00  to  2400  cm  ^ in  steps  of  5 cm  . The  spectral  resolution 
is  5 cm  For  both  species  and  both  parameters,  a linear  extrapolation  of 
log  k or  log  6 with  T is  used  to  obtain  values  at  T > 3000  K. 

The  line  width  band  model  parameter  for  pressure  (Lorentz)  broaden- 
ing is  taken  from  the  NASA  handbook,  but  simplified  to  allow  only  a single 
foreign  gas  broadening  component.  The  form  is 


Y (v,  c,  p,  T) 


= Yo(v)p[b,  c(^)  + c’l/M 


+ a2^^  ■ 


(142) 
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Yq(v)  is  the  broadening  parameter  at  standard  temperature  and  pressure  for 
non-resonant  self-broadening,  is  the  efficiency  (relative  to  Yq)  for  resonant 
self-broadening  and  ct^  is  the  efficiency  for  foreign  gas  broadening.  These 
efficiency  factors  are  assumed  to  be  independent  of  v and  are  aj,  = 6.  53, 

0.2  - for  H2O  in  the  2.  7 -pm  region  and  aj  = 0.  Ill  ^2  ~ 0*  ”78  for  CO2 
in  the  4.  3 -pm  region.  7q(v)  for  the  two  species  is  tabulated  in  Ref.  23  for 
the  same  spectral  grid  and  resolution  as  k and  S 
For  velocity  (Doppler)  broadening, 


(143) 
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where  v is  in  cm  , T in  ® K,  and  oij  is  the  molecular  weight  of  the  active 
species  in  a.  m.  u.  For  H2O,  ^3  = 1®  and  for  CO2,  o^  - ^4. 

For  use  in  the  present  work,  these  param.eters  are  convolved  with  an 
instrument  filter  function  in  order  to  obtain  parameters  appropriate  to  various 
simulation  and  measurem.ert  conditions.  Let  f(v)  be  a general  filter  function 
that  falls  to  zero  at  Vj  and  ('.  j < normalized  such  that 


/ 


f {•.,  ) d ; = I 


Within  the  statistical  band  model  formulation,  the  appropriate  effective  param- 
eters for  the  bandpass  f(,:)  are  defined  by  (see  Appendix  C) 


k(T)  = 


f 


f(, ) k(.,,  T)  d.. 


(144) 


'1 
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Three  H2O  and  two  CO2  filter  bandpasses  are  used  in  the  current  work. 
Two  of  the  H2O  bandpasses  are  the  actual  filters  used  in  the  NERD  emission/ 
absorption  measurements  program  and  are  illustrated  in  Fig.  8.  f(v)  is 
obtained  by  properly  normalizing  these  filter  transmittance  curves.  The 

third  HjO  filter  is  used  for  simulation  work  and  is  a simple  rectangular 

^ -1  - 1 
function  of  width  Av  = 25  cm  centered  on  m = 4000  cm  . The  CO2  filter 

used  in  the  NERD  measurements  program  is  shown  in  Fig.  9.  A second  CO2 

bandpass  used  in  simulation  calculations  is  rectangular  with  width  Av  = 25  cm 

and  centered  on  v = 2272.7  cm 

The  effective  band  model  parameters  k and  6 obtained  from  (144)-(146) 
for  the  three  H2O  filters  are  shown  in  Figs.  10  and  11.  The  effective  param- 
eters for  the  two  CO2  filters  are  shown  in  Figs,  12  and  13.  The  Yq  and  v 
parameters  for  the  five  filters  are  tabulated  in  Table  3. 
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Table  3.  Pressure  and  Doppler  Line  Broadening  Parameters 
for  the  H2O  and  CO2  Filter  Bandpasses 


Filter  Bandpass 

Yq  (cm“Vatm) 

V (cm'S 

NERD  H2O  Narrowband 

0.0739 

3865.  5 

NERD  H2O  Wideband 

0. 0739 

3985.0 

H2O  Rectangular 

0.0674 

4000. 0 

NERD  CO2 

0.0958 

2282.5 

CO2  Rectangular 

0.0941 

2272.7 

VI.  EMISSION/ ABSORPTION  PROFILE  PREDICTIONS 


Program  EMABIC  was  used  to  generate  transverse  radiance  and  absorp- 
tance  profiles  for  twelve  engine  test  cases.  These  cases  were  chosen  to 
represent  a wide  range  of  engine  conditions  (e.  g. , tVirust,  pressure,  nozzle 
contour).  For  all  cases,  the  radial  profiles  of  pressure,  temperature,  and 
concentration  (pTc  profiles)  at  the  nozzle  exit  plane  were  obtained  with 
standard  combustion /nozzle  codes.  The  engine  designations,  t^pe  ot 

exit-plane  calculation,  and  case  numbering  are  listed  in  Table  4.  For  the 
first  seven  cases  (the  1 -dimensional,  constant  oxidizer -to -fuel  ratio  cal- 
culation cases),  the  source  conditions  are  uniform  across  the  exit  plane. 

These  source  conditions  are  tabulated  in  Table  5.  The  last  five  cases 
represent  nozzle  exit  planes  with  varying  degrees  of  nonuniform  source 
conditions.  The  pTc  profiles  for  these  cases  are  shown  in  Figs.  14- IS. 

Transverse  profiles  were  generated  for  both  H2O  at  2.5-pm  and  CO2  at 
4.4-pm  with  the  parameters  for  the  25  cm”^  wide  rectangular  filter  band- 
passes  discussed  in  Section  V.  Calculations  were  performed  for  both  the 
Lorentz  and  Voigt  lineshapes  and  (for  the  nonuniform  source  cases)  for  both 
the  CG  and  DR  nonuniformity  approximations.  Voigt  lineshape  calculations 
were  performed  with  the  RW  approximation.  A spatial  resolution  correspond- 
ing to  N = 50  zones  was  used  in  all  cases.  The  results  for  the  transverse 
radiance  and  absorptance  profiles  are  presented  in  Figs.  19-30. 

For  two  of  the  test  cases  (Cases  3 and  11),  profiles  were  generated 
with  the  modified  NASA  Voigt  line  model  also.  A comparison  of  the  results 
obtained  for  the  two  models  is  shown  in  Figs.  31  and  32.  The  difference 
between  the  models  amounts  to  only  ~ 10%  in  radi,  .ice  and  absorptance. 

The  primary  reason  for  the  generation  of  these  profiles  was  to  obtain 
a familiarity  with  the  magnitudes  and  shapes  of  profiles  that  are  likely  to 
occur  in  test  measurements  on  typical  engines.  Some  of  the  significant 
features  will  be  pointed  cut  in  latter  sections  where  inversions  of  a limited 

These  profiles  were  provided  by  J.  D,  Stewart,  AFRPL. 
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Table  4.  Engine  Test  Cases  for  Emission /Absorption 
Profile  Predictions 


Engine  Designation 

1-D 

Constant  O/F 

2-D 

Constant  O/F 

2-D 

Variable  O/F 

CASE 

CASE 

CASE 

Titan  Il/Stage  1 

1 

Titan  lIIC /Stage  2 

2 

8 

10 

Transtage  Baseline 

a 

11 

ERASE  i.OOO  lb  (6:1  Cone) 

4 

ERASE  500  lb  (25:1  Cone) 

5 

ERASE  500  lb  (25:1  Bell) 

6 

9 

12 

300  lb  MM  Bus 
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Rad 


Radial  Profiles  for  Engine  Case  9.  (a)  Temperature;  (b)  Pressure 

ic)  H?0  and  CO2  concentrations. 


Fig. 


Fig.  18.  Radial  Protiies  tor  JbJngine  Case 
(c)  H2O  and  CO2  concentrations. 


Fig.  ZO.  Transverse  Profiles  for  Engine  Case  2.  (a)  H2O  Radiance;  (b)  H^O 

A.bsorptance;  (c)  CO2  Radiance;  (d)  CO2  absorptance. 


Radianc 


e Cas 


e Cas 


Profiles  for  Engine  Case  9.  (a)  H2O  Radiance;  (b)  H^O 

; (c)  CO2  Radiance;  (d)  CO2  absorptance. 


Radiance;  (b)  H2O 


Fig.  29.  Transverse  Profiles  for  Fngine  Case  11.  (a)  H2O  Radiance;  (b)  H2O 

Absorptance;  (c)  CO2  Radiance;  (d)  CO2  absorptance. 


Fig.  30.  Transverse  Profiles  for  Engine  Case  12.  (a)  H2O  Radiance;  (b)  H^O 

Absorptance;  (c)  CO2  Radiance;  (d)  CO2  absorptance. 


z/R 
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Fig.  31. 


Comparison  of  Transverse  H^O  Profiles  for  Case  3 Generated  with  the 
RW  and  Modified  NASA  Voigt  Line  Models,  (a)  Radiance;  (b)  Absorptance 
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VII.  INVERSION  CHARACTERISTICS 


A.  Introduction 

The  intrinsic  inversion  characteristics  of  the  iterative  Abel  inversion 
algorithm  are  considered  in  this  section.  The  basic  method  involves:  (1)  the 
assumption  of  radial  pTc  profiles;  (2)  generation  of  transverse  emission/ 
absorption  profiles  from  these  radial  profiles;  (3)  an  immediate  inversion  of 
these  transverse  profiles  to  retrieve  the  radial  profiles,  cind  (4)  a comparison 
of  the  retrieved  and  true  radial  profiles.  No  consideration  is  made  here  on 
the  effect  of  random  or  bias  error  in  the  "synthetic  data"  (the  transverse 
profiles  generated  in  step  2)  on  inversion  (except  for  that  introduced  by 
rounding  the  output  of  step  2 to  a finite  number  of  significant  digits  before 
inserting  it  into  step  3).  Thus,  this  analysis  assesses  only  the  convergence 
and  uniqueness  aspects  of  the  inversion  model.  Some  effects  of  random  and 
bias  data  error  are  treated  in  Sections  VIII  and  IX,  respectively. 

Three  test  cases  are  treated  in  the  present  analysis,  lue  first  is  a 
purely  hypothetical  engine  case.  The  second  and  third  are  the  engine  test 
cases  4 and  12  considered  in  the  previous  section. 

B.  Hypothetical  Engine  Test  Case 

Two  cases  involving  hypothetical  radial  profiles  were  considered.  The 

assumed  radial  pTc  profiles  for  the  first  (Test  Case  1)  are  shown  in  Fig.  33. 

The  radial  profiles  for  the  second  (Test  Case  2)  are  identical  in  shape,  but 

the  nozzle  radius  is  increased  an  order  of  magnitude  from  5 to  50  cm  and  the 

concentration  is  increased  from  0.3  to  1.0.  For  both  cases,  the  assumed 

active  species  is  H,0,  and  the  25  cm~^  wide  rectangular  filter  centered  on 
- 1 ^ 

4000  cm  is  used.  Transverse  emission /absorption  profiles  were  generated 
with  a spatial  resolution  grid  corresponding  to  N = 20  zones,  the  Lorentz 
lineshape,  and  the  CG  nonuniformity  approximation,  (All  of  these  calcula- 
tion conditions  were  duplicated  in  inversion).  The  results  are  shown  in  Fig,  34. 

A slightly  different  line -broadening  model  from  that  described  in  Sevition  V 
was  used  in  these  calculations,  but  this  sanie  model  was  used  in  the  subsequent 
inver  sions. 
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Hypothetical  Engine  Test  Cases. 
;k  source). 


From  the  degree  of  centerline  absorptance  for  these  two  cases,  we  can 
describe  Test  Case  1 as  representing  a "relatively  thin"  source  condition 
and  Test  Case  2 as  a "moderately  thick"  source  condition. 

Inversion  of  Test  Case  1.  Figure  35  shows  the  Abel  inversion  results  for 
Test  Case  1.  Since  the  source  is  relatively  thin,  these  results  are  fairly 
close  to  the  "true"  radial  profiles. 

A measure  of  the  improvement  upon  the  simple  Abel  inversion  by 
successive  iterations  is  shown  in  Fig.  36.  Figure  36a  is  a plot  of  the  root- 
mean- square  (rms)  difference  between  successive  iterations  for  T and  c for 
increasing  iteration  number.  If  f is  either  T or  c,  then 


&f 


rms 
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The  curve  for  AT  indicates  convergence  near  i = 6 with  a final  convergence 
residual  of  AT  0.2  K.  The  concentration  curve  indicates  convergence 
near  i 11  with  a residual  value  of  Ac  =»■  b x 10~^.  Both  of  these  residual 
values  are  very  small  compared  to  the  absolute  values  of  the  true  T and  c 
profiles.  It  is  conjectured  that  finite,  rather  than  zero,  values  of  the  resid- 
uals occur  because  only  four  significant  figures  are  retained  from  the 
generated  transverse  profiles  for  input  into  the  inversion  process. 

A more  significant  test  for  convergence  is  illustrated  in  Fig.  36b. 

Here,  the  rms  differences  between  the  input  transverse  profiles  and  trans- 
verse profiles  generated  from  successive  iterations  for  T and  c are  plotted. 
Both  the  curves  for  radiance  (N)  and  transmittance  (t)  show  a convergence  at 
i ^ 5 with  residuals  of  AN  3 x 10"^  and  aT  1.  5 x 10"^. 

The  detailed  difference  between  the  inverted  T and  c profiles  and  the 
true  T and  c profiles  are  shown  for  i = 5 in  Fig.  37a.  Note  that  this  display 
is  not  related  to  Fig.  36a  which  was  a measure  of  difference  between 
successive  iterations  for  T and  c.  The  deviations  of  Fig.  37a  display  a 
random  distribution  of  positive  and  negative  deviations  in0srs:R=5  cm. 
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The  magnitude  of  the  deviations  also  appear  to  be  normally  distributed. 

The  only  conspicuous  behavior  is  a tendency  for  the  maximum  deviations 
to  occur  at  the  end  point  r = R.  This  behavior  can  be  traced  to  the  numerical 
differentiation  algorithm  of  the  Abel  inversion  routine.  Over  most  of  r,  a 
central -difference  formula  is  used;  at  r = R,  a backward  difference  formula 
is  used. 

Figure  37b  is  similar  to  Fig.  37a,  but  shows  the  detailed  difference 
between  the  input  transverse  profiles  and  the  transverse  profiles  generated 
from  the  fifth  iteration  for  T and  c.  A summary  of  Fig.  37  is  given  in  Table 
6.  The  totality  of  these  results  indicate  a unique  convergence  to  the  true 
radial  profiles 

Inversion  of  Test  Case  2.  Simple  Abel  inversions  of  the  transverse  profiles 
of  Fig.  34b  yield  the  T and  c profiles  of  Fig.  38.  Here,  because  the  source 
is  not  even  approximately  thin,  the  Abel  inversion  results  are  not  at  all  close 
to  the  true  radial  profiles.  Convergence  effected  by  successive  iterations  is 
shown  in  Fig.  39.  Figure  39a  is  similar  to  Fig.  36a  for  Test  Case  1 and 
shows  the  convergence  between  successive  radial  profiles.  Figure  39b  is 
similar  to  Fig.  36b  for  Test  Case  1 and  shows  the  convergence  between  the 
input  transverse  profiles  and  profiles  generated  from  successive  c and  T 
profiles.  Both  Figs.  39a  and  39b  indicate  convergence  near  i = 30.  Details 
of  the  radial  and  trjinsverse  profiles  at  i = 30  are  shown  in  Fig.  40.  The 
discussion  of  this  figure  is  the  seune  as  the  discussion  of  Fig.  37  for  Test 
Case  1.  A summary  of  the  results  of  Fig.  40  is  given  in  Table  7.  As  for 
Test  Case  1,  these  results  indicate  a unique  convergence. 

The  only  significant  difference  between  the  Test  Case  1 and  Test  Case 
2 inversions  is  the  number  of  iterations  required  for  "ultimate"  convergence, 
that  is,  convergence  to  the  point  where  the  transverse  profiles  generated 
from  successive  iterations  for  c and  T do  not  become  closer  to  the  input 
transverse  profiles.  For  Test  Case  1,  five  iterations  were  required  while 
for  the  thicker  scarce  of  Test  Case  2,  thirty  iterations  were  required.  Al- 
though this  latter  number  may  appear  excessive,  the  inversion  routines  are 
reasonably  efficient,  and  thirty  iterations  are  completed  in  less  them  1.5 
seconds  of  computing  time  on  the  CDC  7600  computer. 
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Table  6.  Summary  of  Inversion  Characteristics  for  Test  Case  1 


Variable 

Mciximum 

Profile 

Value 

Deviation  Results  at  i = 5 

Mean 

* 

rms 

Maximum 

T 

2000  K 

0.33 

3.5 

4.0 

c 

0.30 

0. 0024 

0.0021 

0.0068 

N 

1 . 5(-4)W/cm^-8r-cm”  ^ 

2.9(-10) 

3.3(.8) 

7.4(-8) 

1 -T 

0,044 

3.3(-7) 

1.7(-5) 

3.4(-5) 

from  true 

except  endpoint  at  r = R 
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Table  7.  Summary  of  Inversion  Characteristics  for  Test  Case  2 


Variable 

Maximum 

Profile 

Value 

Deviation  Results  at  i = 30 

Mean 

rms* 

Maximum 

T 

2000  K 

0.43 

2.4 

4.0 

c 

1.00 

0.0012 

0.0071 

0.026 

N 

2.  0(-3)W/cm^-8r-cm”^ 

1.5{-8) 

3.3(-7) 

7.0(-7) 

1 - "t 

0.45 

-1.6(-6) 

3.2(-5) 

7.4(-5) 

from  true 

except  endpoint  r = R 


-115- 


C.  Engine  Test  Case  4 

Engine  Test  Case  4 represents  a uniform  source  with  R = 4.  178  cm, 
p = 0.790  atm,  T = 1728  K,  c„  _ = 0.373,  and  c_^  = 0.0597.  Transverse 

vj  Cy  Uo 


I 

I 

f 


profiles  generated  from  these  conditions  for  the  Voigt  lineshape  (RW  approx- 
imation) and  N = 50  are  shown  in  Fig.  22.  Inversions  of  these  profiles  were 
performed  with  these  same  conditions.  Abel  inversion  results  for  T and  c 
are  shown  in  Fig.  41.  The  tei'nperature  profile  retrieved  'rom  the  H2O  data 
at  2.5-pm  is  very  close  to  the  true  value  of  1728  K because  the  source  is 
optically  thin  in  this  spectral  region.  The  profile  obtained  from  inversion  of 
the  CO2  profiles  at  4.4-pm  is  substantially  different  from  the  true  value 
because  the  source  is  not  thin  in  this  spectral  region.  Similarly,  the  Abel 
inversion  for  the  H2O  concentration  profile  is  closer  to  the  true  H2O  profile 
than  the  result  for  the  CO2  profile  is  to  its  true  profile. 

The  convergence  of  the  inverted  profiles  in  the  iterative  Abel  scheme  is 
shown  for  H2O  in  Fig.  42.  Figure  42a  displays  the  convergence  of  successive 
radial  profiles  and  Fig.  42b  shows  the  convergence  of  successive  transverse 
profiles  to  the  input  profiles. 

These  results  are  qualitatively  similar  to  those  of  the  hypothetical 
engine  test  case  and  again  confirm  the  convergence  and  uniqueness  of  the 
iterative  inversion  scheme. 

D.  Engine  Test  Case  12 

Engine  Test  Case  12  represents  a highly  nonuniform  radiation  source. 
The  radial  pTc  profiles  are  given  in  Fig.  18  and  the  calculated  transverse 
profiles  in  Fig.  30.  The  transverse  profiles  generated  with  the  Voigt  line- 
shape  and  DR  nonuniformity  approximation  were  used  as  input  to  the  inversion 
(and,  of  course,  inverted  with  these  same  conditions).  Inversion  results  for 
temperature  and  concentration  are  shown  in  Figs,  43  and  44,  respectively. 

Explanations  for  the  results  in  the  region  below  r/R  s 0.  2 and  r/R  > 0.  8 
are  required.  Near  the  edge  of  the  source,  both  the  temperature  and  concen- 
tration profiles  plunge  rapidly  to  low  values  (see  Fig.  18).  Consequently,  the 
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CASE  12 

NO  RANDOM  ERRORS 


(b)  Iterative  Abel  Inversion. 


Fig-  44,  Concentration  Inversions  for  Engine  Test  Case  12,  (a)  H2O;  (b)  CO^ 


resulting  transverse  profiles  of  radiance  and  transmittance  are  nearly  zero 
and  unity,  respectively,  near  the  edge  (above  ~0.  80)(eee  Fig.  30).  Since 
essentially  no  information  is  contained  in  this  region,  no  inversion  for  radial 
profiles  can  be  made.  The  inversions  were  carried  out  by  redefining  the 
source  radius  to  be  0.84  of  the  true  source  radius. 

Below  r/R=»'  0.2,  the  temperature  and  concentration  again  plunge  to 
low  values.  Here,  the  interpretation  of  the  inversion  problem  teikes  several 
forms.  An  important  effect  in  the  present  case  is  that  the  transverse  profiles 
were  generated  with  only  four  significant  digits  in  radiance  and  transmittance. 
Since  the  transmittance  is  so  high,  the  effect  of  using  only  this  number  of 
significant  digits  is  to  produce  a random  error  in  the  input  data  simply 
through  the  process  of  round  off. 

However,  results  obtained  with  the  retention  of  many  more  significant 
digits  displayed  the  same  spurious  results  in  this  region.  An  explanation  of 
the  results  can  be  made  with  the  physical  picture  provided  by  the  "onion-peel” 
method.  As  we  work  in  from  the  outer  edge  and  determine  T and  c for  each 
concentric  zone,  we  reach  a point  where  we  must  start  to  subtract  two  leurge 
radiance,  values  in  order  to  determine  a small  radiance  (because  temperature 
is  low)  caused  by  an  inner  zone.  An  inspection  of  the  emission  function  J (r) 
indicates  that  in  order  to  perform  this  subtraction  accurately,  more  than  ~20 
significant  digits  must  be  retained  throughout  the  inversion.  This  interpreta- 
tion is  consistent  with  the  interpretation  that  inversion  can  be  performed  only 
if  the  transverse  profiles  maintain  a spatial  gradient  sensibly  different  from 
zero.  Inspection  of  case  12  profiles  show  that  the  transverse  profiles  are 
essentially  constant  for  r/R  :S  0.  1.  (Note  that  this  same  interpretation  can  be 
applied  to  r/R  ^ 0.  8).  This  feature  is  not  lo  be  interpreted  as  a failing  of  the 
inversion  method,  but  rather  a confirmatic  n of  the  fact  that  an  inversion  can 
be  made  only  if  information  is  available. 

Outside  of  these  two  regions,  the  application  of  the  iterative  inversion 
method  provides  adequate  inversion  convergence. 
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VIII.  RANDOM  ERROR  ANALYSIS 


A.  Introduction 

The  effects  of  random  experimental  error  on  the  accuracy  of  inversion 
was  investigated  for  the  three  test  cases  considered  in  the  previous  section. 
The  approach  is  the  same  as  used  to  investigate  the  intrinsic  inversion 
characteristics:  transverse  profiles  were  generated  from  known  radial 
profiles,  inverted  to  obtain  radial  profiles,  and  the  results  compared  with 
the  starting  radial  profiles.  Here,  the  added  steps  were  taken  of  super- 
imposing artificial  random  fluctuations  onto  the  transverse  profiles  and 
smoothing  these  profiles  before  performing  the  inversion.  Since  the  inver- 
sion has  been  determined  to  be  convergent  and  unique  in  the  absence  of  random 
fluctuations,  any  significant  discrepancy  between  the  inverted  and  starting 
radial  profiles  can  be  attributed  directly  to  the  random  errors  or  the 
smoothing  routine  used  to  suppress  them  or  both. 


B.  Hypothetical  Engine  Test  Case 


The  moderately  deep  (Test  Case  2)  test  case  considered  in  Section 

VIIA  was  used  for  this  random  error  analysis.  The  assumed  radial  profiles 

are  shown  in  Fig.  33  (except  that  the  nozzle  radius  is  R = 50  cm  and  c^j  q = 1). 

2 


The  transverse  emission /absorption  profiles  are  shown  in  Fig.  34b  and  again 
here  in  Fig.  45.  Random  errors  with  a ±5%  rms  deviation  of  the  peak  value 
were  superimposed  on  both  the  radiance  and  absorptance  profiles  at  the  N = 20 
transverse  zone  positions  in  order  to  simulate  experimental  random  fluctua- 
tions. The  magnitude  of  the  individual  fluctuations  was  obtained  from  tables 
of  random  normal  numbers,  and  different  tables  were  used  for  each  profile 
so  that  the  errors  would  be  uncorrelated  between  the  profiles.  These  treuis- 
verse  "error  profiles"  are  shown  in  Fig.  45  along  with  the  true  transverse 
profiles. 

The  results  of  variable  pre-smoothing  of  the  data  according  to  the 
routine  described  in  Section  HID  are  shown  in  Fig.  46  for  the  transverse 
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of  Hypothetical  Engine  Test  Case  Z,  (a)  Radiance;  (b)  Absorptance. 


radiance  and  absorptance  profiles,  respectively.  Each  figure  is  a plot  of 

the  rms  difference  between  the  smoothed  and  unsmoothed  profile  o versus  the 

smoothing  parameter  Y.  A criterion  for  selecting  Y is  to  choose  that  value 

which  makes  o about  equal  to  the  estimated  rms  fluctuations  in  the  data  e,. 

The  estimated  rms  fluctuations  in  radiance  is  iO  W/cm  -sr-cm  (i.  e. , 

5%  of  peak  value),  and  the  appropriate  value  of  Y is  seen  to  be  ~16.  For 

-2 

absorptance,  the  estimated  rms  fluctuations  is  10  and  the  corresponding 
Y value  is  ~2.  The  smoothed  profiles  resulting  from  these  choices  for  Y 
are  shown  in  Fig.  45  along  with  the  true  and  error  profiles. 

Finally,  the  radial  profiles  resulting  from  an  inversion  of  the  pre- 
smoothed transverse  profiles  are  shown  in  Fig.  47,  The  temperature  profile 
is  quite  respectable  and  displays  only  a 81  K rms  deviation  from  the  "true” 
profile.  The  concentration  profile,  on  the  other  hand,  is  not  as  good  and  dis- 
plays an  rms  difference  from  the  "true"  concentration  profile  of  over  20%. 
These  results  are  based  on  the  use  of  Y = 16  and  2,  respectively,  for  radiance 
and  absorption.  An  attempt  was  made  to  further  smooth  the  concentration 
result  by  increasing  the  Y value  for  absorption,  but  little  was  gained.  Since 
the  inversion  for  concentration  is  not  governed  solely  by  the  absorption 
profile  but  is  coupled  with  the  radiance  profile,  an  increased  smoothing 
of  the  radiance  profile  may  also  be  needed  in  order  to  smooth  the  concentra- 
tion result. 

A second  analysis  of  this  case  was  run  under  somewhat  different  con- 
ditions. The  spatial  grid  resolution  was  reduced  to  1 cm  by  use  of  M = 50 
zones,  and  random  errors  of  ±.10%  rather  than  ±5%  were  superimposed  on 
the  transverse  radiance  and  absorptance  profiles.  Also,  the  example  was 
run  with  the  consistent  use  of  the  DR  rather  than  CG  nonuniformity  approxi- 
mation. The  true,  synthesized  experimental,  and  appropriately  presmoothed 
transverse  radiance  and  absorptance  profiles  are  shown  in  Fig.  48.  The 
comparison  between  exact  and  inverted  radial  profiles  for  temperature  and 
concentration  are  shown  in  Fig.  49.  The  comparisons  apparent  in  these 
latter  tw’o  figures  are  qualitatively  the  same  as  the  comparisons  reported  for 
the  5%  random  error  example.  However,  the  quality  of  the  comparison  is 
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parison  between  True  and  Inverted  Profiles  for  10%  Error 

ysis  of  Hypothetical  Engine  Test  Case  2.  (a)  Temperature  compari 

(b)  Concentration  comparison. 


worse.  This  is  to  be  expected,  of  course,  since  the  input  data  is  at  least 
twice  as  noisy  as  the  ±5%  error  example.  In  the  present  case,  the  rms 
difference  between  the  true  and  inverted  temperature  profiles  is  140® K,  up 
from  81°  K for  the  5%  error  case.  In  addition,  the  temperature  profile 
appears  to  be  over-smoothed  yet  is  based  on  the  same  quantitative  pre- 
smoothing criteria  used  for  the  5%  example.  The  inverted  concentration 
profile,  on  the  other  hand,  displays  the  scime  ~23%  rms  difference  from 
the  true  profile  (c  = 1 ) that  was  obtained  in  the  5%  error  exeimple,  but  here 
the  error  is  all  in  one  direction  (higher)  than  the  true  value  rather  than 
distributed  around  the  true  value. 


Engine  Test  Case  4 


The  true  radial  profiles  for  this  case  are  given  in  Table  5 and  the 
transverse  profiles  in  Fig.  22,  Random  fluctuations  of  ±2%  rms  of  peak 
profile  value  were  superimposed  on  the  radiance  and  transmittance  profiles. 
These  error  profiles  are  shown  in  Fig.  50  for  H2O  and  Fig.  51  for  CO^. 
Although  the  rmi<  error  magnitude  is  only  2%,  the  fact  that  it  is  superimposed 
on  transmittance  corresponds  to  a superposition  of  "#30%  errors  on  the 
absorptance  profile  and  ~6%  on  the  CO^  absorptance  profile.  Presmoothing 
of  these  profiles  gave  the  following  optimum  values  for  Y (that  is,  these 
values  gave  rms  differences  between  the  smoothed  and  unsmoothed  profiles 
equal  to  the  superimposed  errors). 


H2O  radiance  : Y = 0.  58 
H2O  absorptance;  Y = 1.35 


CO2  radiance;  Y = 0. 40 

CO2  absorptance:  Y = 0.  53 


The  CO2  profiles  smoothed  with  these  CO2  values  for  Y are  shown  in  Fig.  51. 
The  smoothed  profiles  for  H2O  shown  in  Fig.  50  were  obtained  with  Y values 
eight  times  larger  than  these  optimum  H2O  values.  The  reason  for  this 
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Fig.  50.  Transverse  H^O  E/A  Profiles  for  Random  Error  Analysis  of  Engine  Test 
Case  4.  {a)  JEwdiance;  (b)  Transmittance. 


increase  in  V was  that  for  smaller  Y,  the  error  curves  were  "till  too  rough 
to  allow  the  inversion  algorithm  to  operate  without  generating  indeterminate 
conditionr  in  the  starting  Abel  inversion. 

The  resulting  inversions  for  T and  c are  shown  in  Fig.  52.  In  all 
cases,  the  results  display  severe  fluctuations  about  the  true  profiles.  The 
magnitude  of  fluctuations  is  so  bad  that,  excepr  for  the  temperature  res\alt 
obtained  from  inverting  the  4. 4-pm  CO2  profiles,  the  results  are  nearly 
meaningless. 

The  reason  for  these  poor  results  is  due  to  the  heretofore  used  smooth- 
ing  criterion  that  the  optimum  Y is  that  for  which  the  rms  difference  between 
the  smoothed  and  vinsmoothed  transverse  profiles  is  equal  to  the  super- 
imposed rms  error.  This  criterion  is  too  severe.  A reanalysis  of  Case  4 
is  mar.d  In  Section  VIIIE. 

D.  Engine  Test  Case  12 

The  radial  pTc  profiles  and  transverse  emission /absorption  profiles 
for  Case  12  are  shown  in  Figs.  18  and  30,  respectively.  The  inversion  of 
this  case  with  superimposed  2%  rms  random  errors  was  carried  out  only 
for  the  4.4-|am  profiles.  In  the  2. 5 -pm  region,  the  transmittance  is  never 
smaller  than  ~0.  98.  Thus,  the  imposition  of  2%  errors  on  transmittance 
converts  to  a 100%  rms  error  in  absorptance.  Since  absorptance  is  the 
basic  data  used  in  inversion,  it  was  felt  that  inversion  of  this  case  would  be 
meaningless.  The  true,  randomized,  and  smoothed  4. 4-pm  transverse 
profiles  are  shown  in  Fig.  53.  The  smoothing  parameters  were  increased 
to  four  times  their  optimum  value  in  order  to  avoid  indeterminate  results. 

Inversion  results  for  the  Abel  and  iterative  Abel  methods  are  shown 
in  Fig.  54.  As  for  the  case  with  no  errors  (see  Section  VIID),  inversion 
results  foj  r/R  below'  0.4  and  above  0.  85  are  meaningless;  otherwise,  tlie 
inversion  display  the  expected  behavior. 

The  iterative  Abel  inversion  displayed  a new  behavior  not  seen  in  any 
previous  inversions.  With  increasing  iteration,  the  rms  difference  between 
successive  radial  profiles  at  first  fluctuates  around  a generally  decreasing 
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value  before  maintaining  a steady  decrease  at  about  i = 10  (Fig,  55),  However, 
the  decrease,  at  least  for  temperature,  continues  only  out  to  »vi  = 20  at  which 
point  an  instability  develops  which  results  in  an  increase  of  with  i. 

The  inversion  results  of  Fig.  54  were  taken  from  the  20th  iteration. 

E.  Considerations  on  Pre smoothing 

The  results  obtained  so  far  indicate  that  random  fluctuations  in  input 
transverse  emission/absorption  profiles  can  produce  serious  error  in  the 
inverted  radial  profiles  even  when  the  random  fluctuations  are  dampened  by 
presmoothing  of  the  transverse  profiles.  The  criterion  applied  so  far  to  the 
degree  of  smoothing  allowed  on  the  transverse  profiles  is  that  the  rms  dif- 
ference a between  the  smoo^^hed  auid  unsmoothed  profiles  should  be  about 
equal  to  the  estimated  experimental  rms  error  € of  the  data  (or,  in  tire  case 
of  the  data  error  simulation  tests,  equal  to  the  rms  value  of  the  artifically 
imposed  fluctuations).  Reconsideration  of  the  smoothing  procedure  suggests 
that  this  criterion  is  too  severe  in  most  case's  and  that  a higher  degree  of 
smoothing  can  be  tolerated  without  introducing  any  significant  error,  that  is, 
without  seriously  "over -smoothing"  the  profiles. 

Consider  the  a versus  y curves  of  Fig.  46.  These  curves  are  typical  for 
most  of  the  smoothing  performed  so  far.  The  shape  of  the  curve  is  linear 
for  small  Y.  As  Y increases,  the  curve  levels  off  somewhat  to  form  a plateau 
region  before  again  increasing  sharply  with  Y.  The  appropriate  values  of  Y 
based  on  estimated  (actually,  in  this  case,  known)  rms  fluctuations  occur  in 
this  plateau  region.  Presumably,  the  shape  of  the  curve  in  the  linear  region 
corresponds  to  increased  smoothing  of  the  random  errors  with  no  distortion 
of  the  underlying  profile,  the  plateau  region  corresponds  to  complete  smooth- 
ing of  the  random  errors  with  no  pro.lle  distortion  (i.  e. , over  smoothing), 
and  the  sharp  increase  of  a with  Y beyond  the  plateau  corresponds  to  an 
increase  smoothing  and  distortion  of  their  underlying  profile.  Although  not 
shovni  in  Fig.  46,  with  increased  smoothing  beyond  the  upswing  at  the  end 
of  the  plateau,  the  curve  again  levels  off  at  a value  near  the  mean  value  of 
radiance  (or  absorptance).  This  condition  reflects  a complete  smoothing 
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(mole  fraction) 


(errors  and  underlying  profile)  to  some  constant  value,  that  is,  to  the  point 
■where  the  curvature  of  the  smoothed  profile  reaches  its  absolute  minimum 
value  of  zero. 

The  usual  condition  that  obtains  in  these  analytic  cases  is  that  the 
value  of  Y that  makes  o equal  to  the  superimposed  rms  error  e occurs  near 
the  left  hand  side  of  the  plateau.  In  application  to  real  data,  e iu  only  an 
estimate  and  such  an  exact  selection  of  Y could  not  be  made..  However,  since 
a does  not  vary  much  with  Y in  the  plateau  region,  its  selection  is  not  critical. 
Thus,  a second  criterion  for  the  selection  of  Y can  be  imposed  by  requiring 
it  to  be  in  the  plateau  region  (if  it  can  be  identified  in  the  o vs  Y plot).  How- 
ever, even  this  is  not  enough.  Consider  the  a versus  y curve  for  the  Case  4 
Engine  Test  shown  in  Fig.  56.  Here,  the  plateau  region  extends  over  an 
enormous  range  from  ~ 1 to  beyond  1000  (the  upswing  is  still  not  e-vident  at 

Y = ICOO).  The  optimum  value  of  Y determined  by  requiring  0 = e = 0,02  is 

Y ~ 1.  But,  0 is  not  much  greater  than  0.02  even  for  Y as  large  as  1000. 
Consequently,  smoothing  of  the  profile  with  Y = 1000  does  not  constitute  a 
significant  over  smooching.  In  Fig.  57  are  shown  the  unsmoothed  transmit- 
tance profile  and  the  smoothed  profile  for  Y ~ 1 and  Y = 1000.  For  either  of 
these  cases,  a ~ 0.02,  yet  the  curve  for  Y = 1000  is  decidedly  smoother. 

The  question  remains,  how  can  two  drastically  different  values  of  Y give 
significantly  different  degrees  of  smoothing  yet  give  the  same  result  for  the 
rms  difference  between  the  smoothed  and  unsmoothed  curves?  The  answer 
lies  partially  in  the  distribution  of  differences  between  the  smoothed  and 
unsmoothed  curves.  When  Y is  taken  from  the  left  end  of  the  plateau,  the 
deviations  making  up  a are  all  about  the  same  size.  When  Y is  taken  from 
the  right  end  of  the  plateau,  the  deviations  are  more  randomly  distributed 
(i.  e. , both  small  and  large  deviations  occur).  Cleourly,  the  latter  case  is 
the  condition  desired  since  this  is  how  the  de'viations  were  superimposed  in 
the  first  place.  The  result  of  this  analysis  is  a new  criterion  for  selecting 
Y;  Y should  be  taken  as  the  rightmost  value  of  the  plateau  region.  This 
ensures  that  a ~ e and  that  the  de-viations  are  randomly  distributed. 
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Fig.  57.  Transverse  Transrr.ittance  Curves  for  Random  Error  Analysis  (2%  rms 
error)  of  Engine  Test  Case  4 


I 
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The  Case  4 inversion  for  H^O  was  rerun  with  increased  smoothing. 

The  critical  variable  here  is  the  degree  of  noise  on  the  transmittance  or 
absorptance  profile.  The  effect  on  increased  smoothing  on  this  profile  has 
already  been  shown  in  Fig.  57.  Inversions  were  made  for  smoothing  param- 
eters of  Y-  = 10,  100,  1000  (the  value  Y = 10  is  near  the  value  used  in  the 
a.  ^ 

inversion  of  Section  VIIB).  The  radiance  profile  was  smoothed  with  corre- 
sponding values  Vj.  = 5,  50,  100. 

The  temperature  inversion  results  for  these  increasing  values  of  Y 
are  shown  in  Fig.  58a.  The  large  spikes  for  the  previc,  -■  inversion  (V^  = ’.0) 

C. 

are  substantially  reduced  for  = 100  and  are  essentially  nonexistent  for 
Y = 1000.  In  the  last  case,  the  maximum  deviation  of  the  inverted  tempera- 
ture  profile  from  the  true  profile  is  only  ~ 200  K.  The  inversion  results  for 
H^O  concentration  are  shown  in  Fig.  58b  where  again  the  enhancement  of 
smoothing  results  in  ^ much  better  inversion,  although  not  as  good  as  for 
temperature.  Here,  even  for  y = 1000,  a 30%  discrepancy  occurs  between 

cL 

the  inverted  and  true  profile.  However,  this  discrepancy  cannot  be  critisized 
since  the  errors  of  the  original  absorptance  curve  are  also  ~ 30%. 
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IX.  BIAS  ERROR  ANALYSIS 


A.  Introduction 

A variety  of  calculations  were  performed  in  order  ^o  assess  the  effects 
of  bias  error  on  the  accuracy  of  radial  pTc  profile  retris  /al.  With  the  excep- 
tion of  assuming  an  infinitesimal  field  of  view  on  data  collection,  bias  error 
in  the  input  data  itself  (e.g.  , errors  due  to  radiometer  calibration  error, 
instrument  zero  line  drift  and  optical  alignment  errors)  were  not  considered 
because  they  are  generally  very  much  smaller  than  the  random  error  of  the 
data  and  essentially  masked  by  the  latter.  (An  exception  to  this  is  presented 
in  the  next  section  where  the  effect  of  correcting  the  radiance  profile  for 
background  radiation  is  considered.  ) The  bias  errors  considered  here  are 
those  that  may  be  introduced  by  the  radiation  and  inversion  models  used  in 
the  analysis,  that  is,  those  introduced  by  the  modelling  assumptions. 

B.  Effect  of  Band  Model  Parameter  Uncertainty 

The  effect  of  the  uncertainty  in  band  model  parameters  on  inversion 
is  considered  in  two  ways.  First,  a simple-minded  analytical  treatment 
is  presented.  Second,  the  results  of  some  actual  inversion  in  which  param- 
eter variations  were  introduced  are  discussed. 

Analytical  Analysis 

In  any  band  model  for  an  array  of  arbitrary  lines  and  a uniform  optical 
path,  the  weak  and  strong  absorption  limits  for  absorptance  are 


a 


ku  weak 
1 strong 


These  limits  simply  display  the  linear  dependence  of  a on  k for  weak  (a  0) 
absorption  and  the  independence  of  a on  any  band  model  parameter  for  strong 
absorption.  Consider  then,  only  the  weak  absorption  case.  The  fractional 
change  in  a doe  to  changes  in  k is  simply 


I 


i 
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and  similarly,  since  N = N a.  tiie  fractional  change  in  radiance  and  N is 
also 

❖ _ ^ 

AN  _ an  _ Ah 

N*  " N " k 


The  change  in  inversion  temperature  to  changes  in  N is  related  through  the 
Planck  radiation  law.  In  the  extremes  of  low  and  high  temperature  we  have, 
respectively 


N 


C^e 


•Cj/T 


T - 0 

X -»  00 


where  and  are  radiation  constants.  From  these  limiting  relations  we 
find 


AZ 

T 


'an*  T ^ Ak  T 
N*  ^2  k ^2 


T - 0 


CO 


For  all  practical  purposes,  T < C2  = (1.439  K cm"^)v  so  that  the  worst  case 
variation  is 


AT 

Om 

Ak 

T 

IT 
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A similar  result  for  concentration  error  can  be  obtained  from  the 
thin  source  relation  a = ku  = kcpL.  It  is  easily  established  that 

Ac 
c 

Thus,  we  conclude  that  an  estimate  of  maximum  fractional  temperature 
and  concentration  inversion  error  is  of  the  order  of  the  fractional  error 
in  the  absorption  band  model  parameter  k. 

Inversion  Analysis 

The  band  model  parameters  used  in  EMABIC  are  constructed  by 
combining  the  high  temperature  NASA  parameters  with  low  temperature 
parameters  derived  from  the  AFGL  atmospheric  absorption  line  data  com- 
pilation (see  Section  V).  These  three  sets  of  parameters  for  H~0  at 

-1  -1  ^ 

V = 4000  cm  with  Av  = 25  cm  are  shown  in  Fig.  59.  The  effect  of  band 

model  parameter  variation  is  demonstrated  by  inverting  transverse 

radiance  and  absorptance  profiles  generated  with  the  combined  parameters 

using  the  three  parameter  sets.  Case  11  of  the  engine  test  cases  was  used 

since  this  case  displayed  the  highest  degree  of  temperature  variation  across 

the  plume.  Aside  from  the  variation  of  parameter  set,  all  the  inversions 

employed  the  Voigt  line  shape,  the  DR  approximation,  and  N = 50.  The 

original  transverje  profiles  were  generated  with  these  same  conditions. 

The  inversions  for  temperature  are  shown  in  Fig.  60a.  The  garbage 
below  r/R  = 0.45  is  known  to  be  an  implicit  failure  of  inversion  due  to  the 
lack  of  information  content  in  the  flat  trai. averse  profiles  and  should  be 
ignored  here  (see  Section  VIID).  Above  r/R  = 0.45,  the  solid  line  repro- 
duces the  true  radial  temperature  profile  as  expected.  The  inversion  using 
the  line-averaged  parameters  also  reproduces  the  true  profile.  This  result 
is  also  expected  since  up  to  the  highest  temperature  of  the  optical  path 
{~  1100  K),  the  line-averaged  and  combined  parameters  are  essentially 
identical  (see  Fig.  59).  The  inversion  using  the  NASA  parameters  displays 
a large  deviation  from  the  true  profile  in  tlie  region  0.45  r/R  0.  70. 


Ak 

k 
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'ig.  60.  Inversions  for  Engine  Test  Case  11  with  Three  Band  Model  Parameter 
Sets,  (a)  Temperature;  (b)  H2O  concentration. 
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The  temperature  in  this  region  is  roughly  750  K.  Reference  to  Fig.  59 
shows  that  this  is  precisely  the  temperature  region  where  the  NASA  para- 
meters differ  significantly  from  the  combined  parameters.  At  750  K,  in 
fact,  the  two  parameter  sets  vary  by  Ak/k  ~ 0.  3,  The  inversion  error  in 
0.  5 S r/R  s 0.  7 is  also  ~0.  3,  thus  supporting  the  simple  analytical  result 
obtained  earlier. 

The  inversion  results  >r  H2O  concentration  are  shown  in  Fig.  60b. 
Again,  the  results  obtained  with  the  combined  and  line-averaged  parameters 
are  comparable  while  the  result  obtained  with  the  NASA  parameter  set  is 
quite  different.  This  difference  is  as  readily  explained  as  for  temperature 
although  the  difference  occurs  in  a slightly  different  r/R  region.  Here,  an 
average  error  of  ~35%  obtains  in  the  0.65  -•  0.85  region.  The  temperature 
is  roughly  900  K and  the  percentage  difference  in  k is  ~25%. 


C.  Effect  of  Pressure  Profile 


f 


As  discussed  in  Section  IIIA,  it  is  not  possible  to  im-srt  the  two  inde- 
pendent E/A  transverse  profiles  co  get  all  three  cf  the  radial  pTc  profiles. 
One  of  the  radial  profiles  has  to  be  assumed.  In  all  of  the  inversions  treated 
so  far,  T and  c have  been  assumed  to  be  the  unknown  radial  profiles.  The 
pressure  profile  has  been  taken  as  its  true  value.  In  a realistic  application, 
p would  not  be  known.  (However,  pressure  is  the  one  variable  that  can  be 
reliably  predicted  a priori  with  standard  combustion/nozzle  codes.  ) A test 
was  made  to  assess  the  importance  of  knowing  p.  This  test  was  carried  out 
on  the  engine  test  Case  4.  The  uniform  source  conditions  are  given  in 
Table  5 and  the  computed  transverse  E/A  profiles  in  Fig.  22.  Both  pro- 
file generation  and  inversion  were  performed  with  the  Voigt  line  shape.  In 
inversion,  the  value  of  pressure  was  varied  ± 20%  from  the  >ralue 
used  in  profile  generation.  The  result  v/as  that  when  p was  20%  too 
low,  the  inverted  result  for  c was  20%  too  high.  Similarly,  when  p 
was  20%  too  high,  c was  20%  too  low.  The  variation  on  p had 
no  effect  on  the  inverted  temperature.  The  important  result  is  that  regard- 
less of  the  uncertainty  in  p (at  least  up  to  ± 20%),  the  partial  pressure 
of  the  active  species,  that  is,  the  product  cp  is  retrieved  with  essentially 
no  error.  These  results  are  summarized  in  Table  8. 
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able  8.  Pressure  Variation  Results  (Case  4) 


D.  Effect  of  Lineshape 

Again  using  Case  4 as  an  example,  the  error  that  could  be  introduced 
by  assuming  a radiation  model  based  on  a pure  Lorentz  profile  was  examined. 
In  this  case,  the  transverse  E/A  profiles  generated  with  the  Voigt  profile 
were  inverted  with  the  Lorentz  profile.  The  resnlts  are  tabulated  in  Table  V. 
As  for  the  effect  of  pressure  profile  variation,  the  variation  of  lineshape  had 
no  effect  on  the  result  for  temperature  retrieval.  The  result  for  H^O  con- 
centration was  4.  2%  higher  than  the  true  value,  and  the  CO2  concentration 
was  2.  7%  higher  than  its  true  value.  Case  4 is  a high  pressure  example,  and 
these  small  inversion  errors  would  be  expected.  For  low  pressure  caacs, 
larger  errors  would  occur. 


E.  Effect  of  Nonuniformity  Approximation 

The  effect  of  nonuniformity  approximation  on  E/A  profile  prediction 
is  evident  from  the  nonuniform  source  engine  test  case  results  (Cases  8-12) 
presented  in  Figs.  26-30.  In  all  cases,  the  DR  approximation  predicts  a 
larger  radiance  and  absorptance  than  the  CG  approximation. 

An  assessment  of  the  error  inherent  in  inversion  using  the  CG  or  DR 
radiation  model  was  made  for  Case  10.  This  case  was  chosen  because  of  iis 
one,  well-isolated  and  strong  nonisothermality  feature  at  r/R=^0.35  (see 
Fig.  16a).  The  inversions  were  performed  on  the  2.  B-fjim  transverse  emission/ 
absorption  profiles  generated  with  the  DR  approximation  and  with  no  super- 
imposed errors.  Inversion  of  these  profiles  with  the  DR  approximation  gave 
temperature  and  concentration  profiles  that  agreed  exactly  with  the  true 
radial  profiles.  The  results  obtained  by  inverting  the  DR  transverse  profiles 
with  the  CG  approximation  are  shown  in  Fig.  61.  The  significant  result  is 
that  the  CG  approximation  failed  to  reproduce  the  temperature  dip  at  r/R  0.  35. 
Agreement  with  the  exact  profile  for  larger  r (r/Z  ja  0.5)  is  very  good  and 
adequate  for  smaller  r (r/Z  ^ 0.  1).  The  error  near  r/R  “■0.35  is  ~400  K. 


-152- 


w. 


Fig.  61.  Effect  of  Nonuniformity  Approximation  on  Inversions  for  Engine  '] 
Case  11.  (a)  Temperature;  (b)  H^O  concentration. 


F.  Effect  of  Finite  Field  of  View 

The  assessment  of  FOV  on  inver^’on  was  carried  out  for  the  Case  4 

nozzle  conditions.  These  conditions  are  representative  of  the  ERASE 

100  lb  (6/1  Cone)  test  and  are  uniform  with  R = 4.  178  cm,  p = 0.790  atm, 

T = 1728  K,  Cjj  Q = 0.373,  and  = 0.0597,  CO^  emission /absorption 

2 2 

profiles  were  generated  for  these  conditions  and  then  convolved  with  rec- 
tangular FOV  slits  of  width  Az/R  = 0.04,  0.08,  and  0.  16.  The  true  and 
convolved  profiles  for  radiance  and  transmittance  with  Az/R  = 0.  16  are 
shown  in  Fig.  62.  The  only  discernable  difference  between  the  true  and 
convolved  profiles  lies  in  a region  Az/R  wide  iiround  z/R  = 1.  These  con- 
volved profiles  were  then  inverted  in  order  to  see  how  the  deduced  temperature 
and  concentration  differed  from  the  true  profiles.  The  results  for  temperature 
and  concentration  are  shown  in  Fig.  63.  The  inversions  were  performed  with 
effective  nozzle  radii  R'  = R + Az/2  rather  than  with  the  true  value  R. 

Even  for  Az/R  as  large  as  0.  16,  no  effect  is  seen  on  the  inversion 
for  temperature  within  the  true  boundary  value  z = R.  Even  over  the  extended 
region  z = R + Az/2,  deviations  from  the  true  profile  occur  only  at  the  end 
point,  and  these  are  not  large.  The  effect  of  FOV  on  concentration  is  quali- 
tatively different  than  on  temperature.  Here,  the  effect  propagates  a distance 
Az/2  from  R into  the  true  source  region.  The  slight  overshoot  at  R - Az/2 
amounts  to  a maximum  veilue  of  only  ~ 10%.  At  z = R,  on  the  other  hand,  the 
concentration  for  all  cases  is  approximately  one -half  the  true  value.  The 
conclusion  for  this  one  case  then,  is  that  the  effects  of  FOV  up  to  16%  of  R 
are  inconsequential  to  temperature  inversion  euid  significant  to  concentration 
only  in  the  region  R - Az  / 2 to  R. 

In  a sense,  this  assessment  for  a uniform  source  is  a worst  case 
condition.  If,  for  example,  the  source  were  nonuniform  with  a temperature/ 
concentration  minimum  near  r = 0,  then  the  transverse  emission/absorption 
profiles  would  be  even  flatter  around  z = 0 than  those  for  a uniform  source. 
While  this  condition  causes  serious  problems  in  the  inversion  itself,  it  does 
mean  that  there  is  less  difference  between  the  true  and  FOV  convolved  profiles. 
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CASE  4 


4.4  urn 
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X.  TRANSTAGE  DATA  INVERSION 


A.  Introduction 


Twelve  test  cases  from  the  NERD  Phase  I measurements  program 


were  selected  for  inversion  analysis.  The  E/A  measurements  were  made 


at  AEDC  on  the  Transtage  engine  operating  with  three  different  injector 
designs  and  three  O/F  ratios.  Eight  cases  were  for  CO^  and  four 
were  for  H2O.  The  case  numbering,  engine  operating  conditions,  and 
spectral  bandpass  are  given  in  Table  10.  The  three  principal  operations 
of  data  preparation,  profile  inversion,  and  random  error  analysis  are 
described  below  in  detail  in  Sections  B,  C,  and  D,  respectively.  Revised 
inversions  for  selected  cases  are  discussed  in  Section  E. 


B.  Data  Preparation  and  Smoothing 

The  first  step  in  data  preparation  consisted  of  plotting  the  raw  tran- 
stage E/A  data  so  that  zero  levels  for  radiance  and  absorptance  and  effec- 
tive source  radii  could  be  selected.  These  data  plots  are  shown  in  Figs. 
64  through  75.  The  dashed  lines  of  these  figures  indicate  the  zero  level 
selected,  and  these  values  are  given  in  Table  11.  In  all  cases,  this  level 
was  determined  by  eye  from  the  signal  level  beyond  the  experimental 
source  radius  of  R^  = 60.  8 cm.  Radiance  data  were  adjusted  by  subtract- 
ing the  zero  level  from  all  data  across  the  profile.  Absorptance  data 
were  adjusted  by  a constant  multiplicative  value  that  forced  the  transmit- 
tance to  unity  outside  the  source  region.  For  both  radiance  and  absorp- 
tance and  for  all  twelve  cases,  the  plots  allowed  a selection  of  R = 62  cm 
as  a reasonable  apparent  source  radius. 

For  each  profile,  data  points  for  which  r < 0 (obtained  experimen- 
tadly  because  the  transverse  scanning  mechanism  could  cross  the  center 
line)  and  r > R = 62  cm  were  discarded.  The  profiles  (each  consisting  of 


These  data  were  supplied  by  ARO,  Inc. , in  card  deck  form  and  through 
appropriate  Air  Force  channels. 
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Table  10.  Transtage  Engine  Operating  Conditions. 


Case  No. 


Injector 


species 


spectral  S'iiter 


Baseline 


ITIPCUMZ 


ITIPCUM3 


1.7 

CO2 

NERD  CO2 

2.0 

2.3 

2.0 

H^O 

i 

1 

NERD  H^O  Nr  r row 

1.7 

^^2 

NERD  CO2 

2.0 

2.3 

2.0 

1 

1 

H^C 

NERD  H^O  Wide 

1.7 

CO2 

NERD  CO2 

2.0 

2.3 

2.0 

«2° 

NERD  H2O  Wide 

Table  11.  Zero  Levels  for  Correction  of  Transtage 
£/A  Data 


Case  No. 

AN(W/ cm^-sr-prn) 

Aor 

1 

0.0050 

0.004 

z 

0.0025 

-0.004 

3 

0.0075 

0.004 

5 

0.005 

0 

6 

0.0075 

0.  012 

7 

0.005 

0.008 

9 

0.005 

0.014 

10 

0.005 

0.  004 

11 

0.005 

-0.004 

4 

0.001 

0.002 

8 

0.008 

0.002 

12 

0.003 

0.001 

i 

j 


t 

J 
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~ 2CC  to  300  points  across  the  transverse  position  and  for  an  unequal  grid) 
v^ere  then  smoothed  using  the  general  smoothing  algorithm  of  Section  IIID. 
The  a %rersus  \ plots  for  this  procedure  are  presented  in  Frgs.  76  through 
79.  The  6.0^  cases  for  common  injector  are  shown  together  and  the  three 
H^O  cases  are  shown  toget}t.r  since  these  groupings  display  common 
smoothing  characteristics. 

Discussions  presented  in  Section  VUI  on  the  analysis  of  artificially 
constructed  noisy  data  broiight  out  two  criteria  for  the  selection  of  an 
appropriate  smoothing  parameter:  (1)  Y should  be  selected  so  that  cr  ~ c 
where  e is  an  estimate  of  experinicntal  rr^s  random  error;  and  (2)  Y 
should  be  selected  from  the  right  hand  side  of  the  plateau  region.  Inspec' 
tion  Of  Figs.  76  through  79  indicate  that  whereas  rr.ost  of  the  absorptance 
smoothing  curves  display  a reasonably  well-defined  plateau  (cases  7 and 
9 bcirg  exceptions),  few  of  the  radiance  smoothing  curves  do  (cases  3,  4, 

7 and  possibly  8 being  exceptions).  Thus,  use  of  the  second  criterion 
alone  would  not  be  adequate  for  determining  the  radiance  smoothing  para- 
meters. However,  used  in  conjunction  with  the  first  criterion,  it  is 
usually  possible  to  identify  a hinc  of  a plateau  in  the  radiance  smoothing 
curves  near  a ~ e and  thus  to  make  a plausible  selection  of  Y.  The  esti- 
mated rms  errors  in  the  raw  d£.ta  are  indicated  in  Figs.  76  through  79  in 
order  to  facilitate  the  identification  of  these  quasi-plateaus.  An  initial 
selection  of  Y values  was  made  on  the  basis  of  these  two  criteria  and 
smoothed  curves  were  generated.  These  were  compared  with  the  raw 
data  and  a subjective  assessment  of  the  fit  was  made.  In  some  cases,  a 
noticeable  under-  or  oversmoothing  was  evident.  These  cases  were  rer\in 
with  appropriately  increased  or  d-'creased  smoothing  parameter.  The 
final  selected  set  of  smoothing  parameters  are  tabulated  in  Table  12  and 
comparisons  between  the  final  smoothed  curves  and  raw  data  are  shown 
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Table  12.  Smoothing  Parameters  for  Transtage 
E/A  Data 


Case  No. 

1 

100 

1000 

2 

100 

1000 

3 

100 

1000 

5 

100 

200 

6 

100 

200 

7 

100 

200 

9 

100 

200 

10 

100 

200 

11 

100 

200 

4 

1000 

10,000 

8 

1000 

10,000 

12 

1000 

10,000 
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in  Figs.  64  through  75.  In  all  cases  where  V was  revised  on  the  basis  of 
the  subjective  judgement  of  goodness  of  fit  of  the  smootned  curve  to  the 
data,  the  revised  y was  still  a reasonable  value  according  to  ttte  initial 
selection  criteria.  For  the  most  pan.,  the  smoothed  curves  fit  ttie  data 
quite  well.  Where  the  fit  is  not  as  good  as  one  might  expect  by  eye,  the 
tendency  is  toward  oversmoothing. 

C.  Inversion  Results 

The  smoothed  curves  of  Figs.  64  through  75  were  defined  on  an  equal 
interval  grid  with  N = 100  and  used  as  the  basic  input  for  inversion.  Inver- 
sions were  carried  -!.t  with  the  Voigt  line  shape  (modified  NASA  model), 
the  DR  nonuniformity  approximation,  and  N = 50  zones.  For  each  case, 
the  band  model  parameter  sets  appropriate  to  the  bandpass  notation  of 
Fig.  10  were  used.  The  pressure  profile  used  for  all  inversions  was 
taken  from  combustion/nozzle  calculation  for  O/F  = 2.  0,  * and  ic  shown  in 
Fig.  80.  Iteration  convergence  was  deemed  complete  when  the  rms  tem- 
perature difference  between  successive  iteration  results  was  < 50  K and 
the  concentration  difference  was  <0.001. 


The  inversion  results  are  shown  in  Figs,  81  through  83  for  and 
in  Fig.  84  for  H^O.  Cases  7 and  10  failed  to  yield  an  inversion  because 
of  indeterminate  results  in  the  initiating  Abel  inversion,  and  no  additional 
smoothing  was  applied  to  the  input  profiles  for  these  cases  in  an  attempt 
to  obtain  an  inversion. 

The  general  quality  of  these  inversion  results  is  good  for  r/R>0. 4. 
Below  this  value,  serious  fluctuations  and  gross  inconsistencies  between 
cases  are  often  evident.  This  effect  is  attributed  to  the  flatness  of  tlie 

Data  supplied  by  Lt.  Zlotkowski,  AFRPL. 
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input  transverse  profiles  for  z/R  < 0.  3 and  has  been  discussed  in  Section 

Vim. 

The  results  for  temperature  profile  above  r/R  “•0.4  are  remarkably 
consistent  from  CO  case  grouping  t > ^rouping  and  with  the  H O cases. 
Although  a slight  structure  can  be  discerned  from  the  figures,  the  profile 
is  nearly  flat  with  value  T “■  1400  ± 200  K.  The  results  for  CO^  concen- 
tration above  r/R  “*0.4  are  also  reasonably  consistent.  Here,  however, 
a definite  structure  is  evident.  The  overall  tendency  is  for  a value  of 
c “•  0.  02  near  r/R  =0.4,  a rise  to  c “*0.  07  near  r/R  = 0.  8,  and  a decrease 
to  c ““O,  01  at  the  source  boundary.  Cases  9 and  11  {Fig.  83b)  show  an 
additional  valley  near  r/R  ““0.  7,  and  cases  5 and  6 hint  at  a valley  near 
the  same  point.  H^O  concentration  for  r/R>  0.4  has  a structure  similar 
to  that  of  CO^,.  With  c ~ 0.4  at  r/R  “•  0.  4,  an  increase  to  c ~ 0.  6 at 
r/R  ~ 0.  7 and  a decrease  to  c ~ 0.  05  at  the  boundary  are  evident. 

Although  these  general  tendencies  above  r/R  0.4  are  clear  enough, 
the  quality  of  the  results  is  judged  not  good  enou':''  to  make  any  sure  state- 
ments on  the  v'ariation  of  tadial  profiles  due  to  engine  O/F  or  injector 
change.  Whether  or  not  the  quality  of  these  results  is  good  enough  to 
evaluate  the  adequacy  of  combustion/nozzle  codes  is  the  suuject  of  a sep- 
arate  AFRPL  in-house  project. 

D.  Random  Error  Analysis 

The  random  error  analysis  of  these  Transtage  inversions  was  car- 
ried out  as  in  Section  VIII  by  assuming  radial  pTc  profiles,  generating 
transverse  E/A  profiles,  adding  artificial  random  noise,  smoothing, 
inverting,  and  comparing  the  retrieved  radial  profiles  with  the  assumed 
profiles. 

The  AFRPL  point  of  contact  for  this  effort  is  Dr.  Dwayne  McCay. 
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The  assumed  radial  profiles  for  temperature  and  concentrations  are 
shown  in  Fig.  b5.  (The  pressure  profile  from  Fig.  80  is  retained).  Above 
r/R  •“  0.4,  the  temperature  profile  is  the  mean  value  of  the  actual  Tran- 
stage  inversions  for  both  and  1^2®*  structure  of  the  profile  below 

r/R  “•  0.  4 was  selected  to  be  consistent  with  predictions  of  TDK  nozzle 
codes,  and  the  magnitude  was  selected  to  match  the  profile  above  r/R  “*  0.4. 
The  concentration  profiles  were  constructed  in  a similar  manner,  but  the 
added  constraint  that  range  3 — ►S  (the  range 

of  values  for  the  IZ  engine  test  cases  of  Section  VI)  was  imposed  below 
r/R  0.4.  The  transverse  E/A  profiles  generated  from  these  radial  pro- 
files are  shov-n  in  Fig.  86  and  are  seen  to  be  quite  reasonable  simulations 
of  the  actual  Transtage  data  of  Figs.  64  through  75. 

Tables  of  201  random  normal  numbers  were  generated  and  used  to 
construct  simulated  noise  profiles  from  the  smooth  transverse  profiles. 
Different  tables  were  used  for  radiance  and  absorptance  in  order  to  elim- 
inate any  chance  of  correlated  noise  between  the  profiles.  The  rms  mag- 
nitude of  fluctuation  was  set  at  a fixed  percentage  of  peak  profile  value. 

For  radiance,  absorptance,  H^O  radiance,  and  1^2®  absorptance, 

respectively,  the  percentages  were  3%,  5%,  13%,  and  20%.  These  values 
are  characteristic  of  the  actual  Transtage  data.  These  profiles  were  then 
smoothed  according  to  the  standard  procedures  outlined  earlier.  Both  the 
noisy  data  simulation  and  smoothed  profiles  are  shown  in  Fig.  86.  The 
smoothed  profiles  were  then  inverted  in  order  to  retrieve  the  radial  pTc 
profiles.  The  process  of  noise  addition,  smoothing,  and  inversion  was 
carried  out  five  times  for  each  species,  each  time  with  different  tables  of 
random  normal  numbers.  The  inversion  results  are  shown  in  Fig.  87. 

Qualitatively,  these  results  confirm  the  earlier  suspicion  that  the 
inversion  results  become  less  trustworthy  for  r/R  < 0.4.  All  of  the  evior 
analysis  temperature  and  concentration  inversions  display  an  increasing 
scatter  below  this  level.  The  resxilts  display  the  most  scatter  as 
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would  be  expected  from  the  greater  magnitude  percentage  error  super- 
imposed on  the  synthetic  transverse  data  profiles.  The  rms  scatter  of 
the  inversion  curves  are  shown  in  Fig.  88  at  r/R  = 0,  0.  2,  0.  4,  0.  6,  0.  8 
and  1.  C.  Again,  a sharp  rise  below  r/R  “•  0.4  is  evident. 


Quantitatively,  the  results  do  not  display  as  large  a magnitude  of 
fluctuations  as  is  suggested  by  the  actual  Transtage  inversion  results  of 
Fig.  81.  One  obvious  explanation  is  that  the  error  analysis  treats  ideal- 
ized data  that  does  not  contain  the  subtle  foibles  and  idiosyncrasies  of  real 
data.  Another  possibility  is  that  the  enhanced  scatter  of  the  Transtage 
inversions  is  indeed  due  to  the  variation  of  O/F  and  injector  from  case  to 
case.  The  third,  and  most  likely  possibility  is  that  only  random  error 
has  been  addressed  in  either  the  Transtage  inversions  or  error  analysis. 
Inspection  of  Fig.  81  (particularly  the  three  cases  and  the  CO^  cases 
9,  10  and  11)  reveals  an  oscillatory  variation  of  wavelength  \ ~ 5-10  cm. 
This  noise  has  been  identified  by  ARO  to  be  a deterministic  noise  which 
should  be  removed  (by  a Fourier  analysis)  from  the  data  before  inversion. 
The  fact  that  it  was  not  removed  in  the  present  analysis  nor  simvilated  in 
the  error  analysis  is  most  likely  the  major  cause  of  the  apparent  discrep- 
ancy between  the  observed  and  predicted  inversion  scatter. 


At  the  same  time,  this  argument  does  not  invalidate  the  error 
analysis  results  concerning  random  error.  That  is,  in  the  inversion  of 
a real  Transtage  data  case,  the  rms  error  curves  cf  Fig.  88  provide  a 
reasonable  estimate  of  the  inversion  error  due  to  random  experimental 
error  in  the  measured  transverse  profiles. 


E,  Revised  Inversions 

Three  sets  of  Transtage  data  were  selected  for  further  analysis. 

Two  of  these  cases  are  repeats  of  Cases  2 and  8 considered  above,  but 
with  the  deterministic  noise  removed.  This  noise  has  been  identified  as 
due  to  a "bouncing"  of  the  scanning  arm  that  carries  the  photometer  and 
emission  source  apparatus  across  the  plume.  This  noise  was  removed  by 
a subtraction  of  the  peak  in  the  Fourier  transform  of  the  transverse  pro- 
files at  the  arm-bouncing  frequency.  In  addition,  corrections  to  accoixnt 
for  background  radiation  were  made.  This  background  has  been  attributed 
to  scattering  of  radiation  from  the  engine  nozzle  which  was  inadvertently 
in  the  sensor  field  of  view.  An  empirical  correction  in  time  and  position 
was  applied.  The  third  case  (designated  as  Case  13)  is  a CO^  scan  made 
duru  o the  Phase  III  portion  of  the  AEDC  NERD  measurements  program. 
All  appropriate  corrections  to  these  data  were  also  made. 

The  raw  transverse  scans  (after  corrections)  for  these  three  cases 
are  shown  in  Figs.  89-91  and  the  variable  smoothing  curves  in  Figs.  92- 
94.,  A listing  of  the  smoothing  parameters  selected  is  given  in  Table  13, 
and  the  final  smoothed  curves  shown  in  Figs.  89-91.  All  of  the  original 
Transtage  inversions  employed  a radius  R = 62  cm.  This  constant  value 
was  not  appropriate  for  the  revised  cases.  Here,  the  selections  R = 67, 

65  and  65  cm  were  made  for  Cases  Z,  8 and  13,  respectively.  The  pres- 
sure curve  of  Fig.  80  was  extrapolated  smoothly  to  0.  0337  atm  at  67  cm 
for  these  inversions. 
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Table  13.  Smoothing  Parameters  for  Revised  Transtage  E/A  Data 


Case  No. 

Ya 

2 (revised) 

100 

500 

8 (revised) 

500 

1000 

13 

10 

100 
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The  inversion  results  for  temperature  and  concentration  are  shown 
in  Figs.  95-97.  The  results  for  Cases  2 and  8 are  not  substantially  dif- 
ferent from  the  results  obtained  in  the  previous  section  (compare  Fig.  95a 
with  81a,  95b  with  81b,  96a  with  84a,  and  96b  with  84b;.  The  inversion 
results  for  the  Phase  III  CO^  case  (Case  13,  Fig.  97)  is  also  qualitatively 
similar  to  those  for  the  previous  CO^  cases  and  the  revised  Case  2.  In 
this  case,  however,  the  results  for  r/R  0.  16  are  meaningless  due  to  an 
inversion  instability.  This  instability  is  undoubtedly  due  to  the  valley  in 
the  region  z < 15  cm  evident  in  both  the  transverse  radiance  and  absorpt- 
ance  profiles  (Fig.  91). 


I 


Fig.  64.  Transverse  E/ A Profiles  for  Transtage  Case  1.  (a)  Raw  and 

smoothed  tadiance  data;  (b)  Raw  and  smoothed  absorptance  data. 
Dashed  lines  indicate  zero  levels. 
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Fig.  65.  Transverse  E/A  Profiles  for  Transtage  Case  2.  (a)  Raw  and 

smoothed  radiance  data;  (b)  Raw  and  smoothed  absorptance  data. 
Dashed  lines  indicate  zero  levels. 
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Fig.  69.  Transverse  E/A  Profiles  for  T 
smoothed  radiance  data;  (b)  Ra\ 
Dashed  lines  indicate  zero  leve 
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\ Fig.  71.  Transverse  E/A  Profiles  for  Transtage  Case  8.  (a)  Raw  ar  1 

I smoothed  radiance  data;  (b)  Raw  and  snioothed  ahsorptance  data 

f Dashed  lines  indicate  zero  levels. 
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Fig.  73.  Transverse  E/A  Profiles  for  Transtage  Case  10.  (a)  Raw  and 

smoothed  radiance  data;  (b)  Raw  and  smoothed  absorptance  data. 
Dashed  lines  indicate  zero  levels. 
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Fig.  75.  Transverse  E/A  Profiles  for  Transtage  Case  12.  (a)  Raw  and 
smoothed  radiance  data;  (b)  Raw  and  smoothed  absorptance  data 
Dashed  lines  indicate  zero  levels. 
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Smoothing  Curves  for  Transtage  Cases  5,  6 and  7.  (a)  Radiance 

(b)  Absorptance. 


Fig.  79. 


Fig.  88 


Fig.  93, 
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Fig.  97.  Inversion  ! 
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XI.  SUMMARY  AND  CONCLUSIONS 

The  overall  objective  of  this  work  is  to  clarify  the  problem  of 
whether  or  not  multi -position  E/A  inversion  can  give  pTc  profiles  with 
sufficient  accuracy  for  propulsion  applications  (e.  g. , evaluation  of  com- 
bustion/nozzle code  prediction  capabilities  and  assessment  of  rocket  engine 
performance).  This  objective  is  both  general  enough  and  subjective  enough 
that  it  can  only  be  answered  with  qualification.  The  overall  impression 
obtained  in  this  work  is  that  such  inversions  can  at  least  provide  a valuable 
guide  as  to  the  validity  of  code  predictions  ii  sufficient  care  is  taken  to 
obtain  good  E/A  data  and  U the  application  is  to  certain  types  of  test  engines. 
The  discussions  presented  in  this  section  attempt  to  quantify  this  conclusion. 

The  accuracy  of  E/A  inversion  depends  on:  (1)  the  accuracy  of  the 
mathematical  and  inversion  models;  and  (2)  the  quality  of  the  input  E/A 
data.  The  question  of  data  quality  is  taken  up  later.  Here  we  consider 
the  validity  of  the  models  comprising  the  overall  inversion  scheme,  these 
being  the  band  radiation  model,  the  smoothing  model  and  the  actual  inver- 
sion algorithm.  The  last  of  these  needs  little  discussion.  The  results  of 
Section  VII  clearly  display  the  iinique  and  correct  convergence  character- 
istics of  the  Iterative  Abel  scheme  to  arbitrary  accuracy. 

The  accuracy  of  the  band  radiation  model  depends  primarily  on  the 
accuracy  of  the  fundamental  band  model  parameters,  the  accuracy  with 
which  optical  path  nonuniforinities  are  handled  and  the  accuracy  with 
which  Doppler  broadening  effects  are  modeled.  The  model  used  in  the 
present  work  incorporates  the  most  recent  developments  of  band  model 
theory  and  is  capable  of  predicting  E/A  profiles  for  the  NERD  CO2  and 
H^O  bandpasses  to  an  overall  accuracy  of  15%  or  better  for  most  all  real- 
istic plume  conditions.  In  inversion,  errors  introduced  by  the  radiation 
model  are  of  the  nature  of  bias  errors  and,  to  a first  approximation,  can 
introduce  errors  in  the  inverted  temperature  and  concentration  profiles  of 
the  same  order.  Specific  analyses  for  variations  in  bandmodel  parameters, 
line  profile  and  nonuniformity  approximation  are  presented  in  Sections  DCB, 
DCD,  and  IXE,  respectively. 
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Any  number  of  procedures  can  be  used  for  data  smoothing.  The 
method  employed  here  is  the  rational  approach  of  pre smoothing  the  input 
profiles  by  fitting  the  data  with  the  curve  having  the  smallest  possible  mean 
square  second  derivative  consistent  with  the  constraint  that  the  rms  dif- 
ference between  the  smoothed  and  unsmoothed  profiles  a should  not  be 
much  more  than  the  estimated  rms  fluctuations  of  the  data  e , When 
applied  to  synthetic  data  (that  is,  data  that  has  been  generated  artificially 
by  superimposing  random  fluctuations  on  smooth  profiles),  the  smoothing 
procedure  is  highly  objective;  optimum  smoothing  can  be  effected  without 
recourse  to  subjective  judgment  solely  by  selecting  the  smoothing  para- 
meter y in  the  plateau  region  of  the  a versus  y curve.  When  applied  to  real 
experimental  data,  the  procedure  is  only  quasi-objective.  Often,  the  o 
versus  y curve  fails  to  display  a well-defined  plateau,  and  some  subjective 
judgment  is  often  required  in  order  to  effect  a good  smoothing.  The  error 
introduced  by  smoothing  is  difficult  tc  judge,  and  probably  not  relevant  in 
any  event  since  the  principal  effect  of  smoothing  is  to  suppress  random 
errois  so  that  an  inversion  can  be  performed.  On  the  otlier  hand,  it  is 
evident  that  some  error  must  be  introduced  since  when  the  best  overall 
smoothed  curve  is  determined,  the  fit  at  local  points  on  the  profile  is  not 
as  good  as  what  one  might  think  it  should  be  by  inspection.  An  example  of 
this  phenomenon  is  given  in  Fig.  73b  where  the  best  overall  smoothed  curve 
fails  to  follow  the  peak  structure  at  z 50  cm.  To  the  extent  that  the 
smoothing  of  the  Transtage  data  (Figs.  67-75)  is  typical,  the  smoothing 
procedure  may  introduce  local  errors  of  the  order  of  10%.  A plausible 
solution  to  this  problem  is  an  application  of  the  smoothing  procedure  used 
here  in  a piece-wise  manner  across  the  transverse  profile. 

The  most  detrimental  aspect  of  E/A  inversions  is  the  propagation  of 
random  errors.  Even  the  slightest  of  random  fluctuations  on  the  transverse 
profiles  can  be  amplified  in  inversion  to  such  an  extent  that  the  retrieved 
radial  profiles  are  meaningless.  A smoothing  of  the  input  data  is  almost 
always  required  in  order  to  obtain  realistic  results.  Based  on  the  inver- 
sions considered  in  this  work,  the  following  crude  estimate  of  final  inver- 
sion accuracy  has  been  found;  Let  the  percent  rins  values  of  fluctuation  in 


the  emission  and  absorption  profiles  be  and  e^,  respectively.  Then, 
after  an  optimum  smoothing  to  damp  out  these  fluctuations,  the  inversion 
process  will  reintroduce  errors  so  that  the  temperature  and  concentration 
profiles  will  have  percent  errors  of  the  order  of  and  e^,  respectively. 

There  are  certain  qualifications  to  this  rule.  Most  notably,  the 
error  near  the  center  of  plumes  will  be  much  larger  than  predicted  if  there 
is  a pTc  valley  at  r = 0.  The  result  of  this  valley  is  relatively  flat  trans- 
verse profiles  around  z = 0,  and,  in  extreme  cases,  an  actual  depression 
around  z = 0.  For  most  of  the  inversions  of  this  work,  this  condition 
obtained. 

Aside  from  this  very  crude  error  estimate  and  the  qualification 
mentioned  above,  very  little  can  be  said  about  the  final  error  of  retrieved 
radial  profiles  in  general.  Specific  engine  types  must  be  looked  at  in 
detail.  Some  general  comments  can  be  made  about  the  conditions  that  will 
yield  the  best  inversions.  (1)  The  absorptance  must  not  be  so  small  that 
the  random  experimental  noise  completely  hides  the  signal.  (2)  At  the 
other  extreme,  the  source  cannot  be  so  optically  deep  that  the  absorptance 
is  near  unity.  Then,  the  percent  error  of  the  signal  is  small,  but  there 
is  no  information  in  the  signal.  (3)  The  transverse  E/A  profiles  should  be 
monitonically  decreasing  functions  with  z.  This  condition  ensures  that 
there  is  not  a depression  in  the  pTc  profiles  at  r = 0.  In  thfe  case  of  rocket 
engines  employing  conical  nozzles,  this  is  indeed  the  case  so  long  as  uni- 
form mixture  ratio  conditions  prevail. 

For  the  specific  case  of  the  Transtage  inversions,  more  quantitative 
statements  can  be  made  concerning  the  error  of  inversion.  Consider  the 
revised  Case  2 and  Case  8 inversions  of  Section  XE.  Inspection  of  the  raw 
data  curves  of  Figs.  89  and  90  (and  Fig.  71b  for  H2O  absorptance)  gives 
the  following  rough  measure  of  the  percent  (of  maximum  profile  value)  rms 
fluctuations  in  the  data. 

CO^  radiance  1% 

CO^  absorptance  4% 

H2O  radiance  3% 

H2O  absorptance  20% 
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The  inversions  for  these  two  cases  yield  the  following  mean  radial  values 
for  temperature  and  concentration  (Figs.  95  and  96) 

T ~ 1300  K 
c„^  ~ 0.  03 


while  the  random  error  analysis  of  Section  XD  yields  the  following  errors 
due  to  random  fluct’::<itions  (Fig.  88)  at  r/R  ~ 0.  4 

AT^-,  ~20  K ~ 1.5%  of  1300  K 

CU2 

ATj^  O ~ 40  K ~ 3%  of  1 300  K 

ACco  ~ • 001  ~ 3%  of  0.  03 

Ac„  ^ ~ 0.  05  ~ 13%  of  0.4 
H^O 

These  results  indicate  the  general  usefulness  of  the  rule  that  percent  tem- 
perature error  is  roughly  equal  to  the  percent  radiance  error  and  that  per- 
cent concentration  error  is  roughly  equal  to  percent  absorptance  error 
(actually,  the  correspondence  is  surprisingly  good  for  this  case).  The 
percent  errors  are  smaller  for  r/R  >0.4  and  larger  for  r/R  < 0.4.  The 
worst  case  occurs  at  r = 0 for  which  the  inversion  errors  are  approxi- 
mately (from  Figs.  88,  95  and  96) 


However,  these  errors  drop  quite  sharply  for  r > 0. 
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In  addition  to  these  errors  due  to  random  fluctuations,  the  errors 
due  to  radiation  modeling  and  smoothing  must  be  added.  The  former  is 
approximately  15%  at  worst,  and  the  latter  is  estimated  from  the  goodness 
of  fit  illustrated  in  Figs.  89  and  90  to  be  no  worse  than  ~ 5%. 

The  inversion  results  for  the  revised  versions  of  Cases  2 and  8 are 
repeated  in  Fig.  98  with  the  one  standard  deviation  error  bounds  suggested 
by  this  error  analysis. 
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Fig.  98.  Final  Inversion  and  Error  Results  for  Transtage  Cases  2 and  8. 

(a)  Temperature  retrieved  from  CO2  data;  (b)  Temperature  retrieved 
from  H2O  data;  (c)  CO2  concentration;  (d)  H2O  concentration. 
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APPENDIX  A 

OPTIMUM  SPECTRAL  BANDPASSES 

A.  Selection  of  v and  Av  for  H^O  Measurements 

Two  important  criteria  on  the  selection  of  wavelength  and  bandpass 
for  multiposition  E/A  inversion  are  to  make  sure  that  only  one  active 
species  is  being  observed  and  to  provide  an  optimum  signal  strength. 
Consider  first  the  2.  7 region.  The  NASA  absorption  band  model 
parameter  k for  CO^  and  H^O  at  T = 1500  K are  shown  in  Fig.  A-1. 

This  parameter  is  a measure  of  radiance  and  abosorptance  signal  strength 
for  thin  sources.  Also  indicated  on  the  figure  are  the  NERD  wide-  and 
narrow-bandpasses  and  the  25  cm"^  bandpass  used  in  the  12  engine  test  case 
inversions.  For  H2O  detection,  the  two  criteria  of  nonoverlap  and  signal 
strength  leave  little  option  on  selection  of  v and  Av.  Two  regions  satisfy 
nonoverlapping:  v > 3750  cm  ^ and  v g 3200  cm”^.  The  former  is  more 
sure  since  it  lies  beyond  the  CO^  band  head. 

Typically,  H^O  absorption  for  plumes  of  interest  is  weak.  Therefore, 
we  want  to  choose  a region  of  maximum  k.  Thus,  in  order  to  ensure  non- 
overlap with  CO^  and  to  obtain  maxim.um  k,  consider  only  the  region  above 
3750  cm  The  optimum  v is  the  peak  near  3925  cm  The  position  of 
this  peak  changes  with  temperature  as  sho\Mi  in  Fig.  A-2.  Thus,  for  a 
general  temperature  variation  within  the  plume  source,  an  optimum  v 
cannot  be  selected  a priori.  Some  estimate  of  temperature  gradiant  is 
required  in  order  to  calculate  the  effective  optimum  v,  or,  more 
reasonably,  this  position  can  be  determined  experimentally.  But  even 
here,  the  optimum  value  could  change  with  transverse  scan  position  because 
the  temperature  gradient  along  the  line  of  sight  will  be  different  for  each 
trauisverse  position.  For  the  transtage  inversions,  the  plume  temperature 
profile  at  the  nozzle  is  nearly  constant  with  T ~ 1400  K.  Thus,  from  Fig. 
A-2,  an  optimum  v of  3925  cm"^  is  indicated. 
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Fig.  A-1.  Spectral  Variation  of  Absorption  Band  Model  Paa 
for  H2O  and  CO2  near  2.  7 M-m 


The  selection  of  an  optimum  bandpass  around  v must  involve  a 
trade-off  on  radiance  or  absorptance  optimization.  In  order  to  maximize 
radiance,  a bandpass  that  covers  as  much  of  the  peak  as  possible  should 
be  used.  That  is,  the  optimum  bandpass  for  radiance  is  3800  to  ~ 4200  cm 
In  order  to  maximize  absorptance,  we  want  to  maximize  the  mean  value  of 
k in  the  bandpass.  That  is,  we  should  choose  a narrow  bandpass  ( ~ 20  to 
30  cm  centered  on  3925  cm  The  overall  optimum  bandpass  must  be 
somewhere  between  these  two  extremes.  The  wide  bandpass  H2O  filter 
presently  used  in  the  NERD  measurem  aits  is  a very  good  optimum  filter 

for  radiance  measurements;  it  essentially  covers  the  entire  H,0  spectral 

-1  ^ 
feature  around  3925  cm  and  is  well  isolated  from  the  COy  band.  The 

“2 

presently  used  narrow  bandpass  filter  centered  at  3866  cm  would  be  a 
better  (by  about  50%)  absorptance  measurement  filter  if  moved  to  3925  cm”  . 

The  selection  of  an  optimum  filter  for  sin-ultaneous  radiance  and 
absorptance  measurements  would  require  a detailed  trade-off  study  of 
whether  noisy  absorptance  or  noisy  radiance  measurements  is  more  critical 
in  inversion.  If  absorptance  error  is  the  most  detrimental,  the  optimum 
filter  should  be  weighted  toward  narrow.  If  radiance  noise  is  most  detri- 
mental, the  weighting  should  be  toward  wide.  Similarly,  if  more  importance 
is  placed  on  a good  inversion  for  temperature,  the  band  should  be  wide.  If 
more  importance  is  placed  on  concentration  inversion,  the  banc  should  be 
narrow.  (This  weighting  assumes  that  the  system  noise  is  independent  of 
bandpass, ) 

The  evidence  obtained  so  far  in  the  transtage  measurements  is  that 
absorptance  data  is  more  noisy  than  radiance  data,  and  consequently,  that 
concentration  inversions  are  mere  uncertain  than  temperature  inversions. 

These  results  suggest  that  future  H^O  measurements  be  made  with  a narrow 

-1  ^ -1 
band  I-'  30  cm  wide)  filter  centered  on  3825  cm 


B.  Selection  of  v and  Av  for  CO^  Measuren.ents 

The  complete  spectral  embedment  of  the  2.  T-^m  band  in  the 
H^O  band  precludes  CO^  measurements  in  this  region.  The  logical  region 
is  the  4.3-um  fundamental  CO2  band.  Here,  however,  overlapping  with 
the  CO  band  at  4.  7-n,m  must  be  considered.  The  NASA  absorption  band 
model  parameter  k for  CO^  and  CO  at  T = 1800  K are  shown  in  Fig.  A -3. 

The  NERD  and  test  inversion  bandpasses  are  also  indica.ted.  Although 
both  of  these  bandpasses  overlap  into  the  CO  band,  the  effect  is  probably 
not  important.  Ac..ording  to  code  predictions,  the  ratio  of  CO  to  CO2 
concentration  at  the  nozzle  exit  plane  is  not  likely  to  exceed  ~ 0.  4 for 
storable  missiles  (e.g.  , Titan  II).  The  transtage  predictions  show  — 0.4 
also.  In  this  case,  for  the  NERD  bandpass,  the  ratio  of  CO  signal  to  CO^ 
signal  is  only  ~ 0.  006.  Even  if  the  CO/  CO2  ratio  were  raised  to  2.  5 
(which  might  obtain  for  an  Atlas),  the  signal  ratio  is  only  0.04.  Neverthe- 
less, since  the  CO-  spectrum  is  reasonably  flat,  a move  to  a safer  region 
-1  ^ 

beyond  2300  cm  might  be  prudent. 

In  complete  distinction  from  the  condition  for  water,  the  measurements 
for  CO2  sometimes  show  too  much  signal.  That  is,  the  source  is  optically 
deep,  the  radiance  approaches  the  blackbody  value,  and  the  absorptance  is 
approaches  unity.  As  these  limits  are  approached,  the  information  content 
of  the  data  is  lost.  For  such  cases,  a further  increase  of  v out  to  the  steep 
slope  of  the  band  head  may  be  desirable.  However,  although  k,  and  thus 
radiance  and  absorptance  will  be  less,  an  accurate  determination  of  k in  the 
bandpass  may  be  difficult  because  of  the  steepness  of  the  slope. 

Aside  from  the  possibility  of  trying  to  reduce  the  signal  strength,  the 
selection  of  v and  Av  for  CO,  measurements  is  straightforward.  Assuming 
a rectangular  bandpass  for  simplicity,  it  should  extend  from  2300  to  2375  cm 


-219- 


k (cm  /atm) 


i/(cm’^) 


Fig.  A-3.  Spectral  Variation  of  Absorption  Band  Model  Parameter 
ior  CO2  and  CO  near  4.3  pm 
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APPENDIX  B 

PHILLIPS-TWOMEY  INVERSION  METHOD 

The  method  of  inversion  developed  by  Phillips  (B-l)  and  Twomey 
(B-2)  has  been  adapted  to  the  present  emission /absorption  inversion 
problem.  The  method  has  been  used  extensively  in  the  inversion  of 
atmospheric  sounding  data  to  get  atmospheric  altitude  profiles  of  temper- 
ature and  species  concentration.  A fundamental  feature  of  the  method  is 
the  explicit  account  of  random  errors  in  the  input  data.  Unlike  the  itera- 
tive Abel  inversion  method  which  employs  presmoothing  of  the  input  data, 
however,  account  of  random  errors  is  effected  by  smoothing  the  output 
profiles.  This  smoothing  is  built  into  the  inversion  algorithm  and,  in 
effect,  smooths  the  output  profiles  before  they  are  even  generated.  In  its 
simplest  form,  the  solution  for  the  function  f(r)  defined  by  the  integral 
equation 


\ K(2,r)  f(r)  dr  - g(z) 


(B-l) 


or,  more  accurately,  its  matrix  equivalent 


X 

Ew.  k..  f.  = g. 

i j i i ® j 


(B-2) 


f = (A*  A + yH)'^A*  g) 


(B-3) 


where  f and  g are  column  matrices  with  elements  f.  (i  = 1,  2,  . . . , N + 1) 
and  gj  0 = 2,  ...,  N + 1),  respectively,  A is  a coefficient  matrix  com- 
posed of  elements  a^.^  = w.^  are  the  quadrature  weights  in  the  summa- 

tion approximation  (B-2)  to  the  integral  equation  (B-l),  (see  Section  IIIE 
of  main  text).  A*  is  the  transpose  of  A,  H is  a smoothing  matrix,  and  y is 
a scalar  smoothing  parameter.  The  smoothing  matrix  H is 
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and  is  appropriate  to  the  smoothing  criterion  that  the  curvature  function 

N+1 


be  a minimum  subject  to  certain  restraints.  The  constraint  is  that  the 
(C) 

profiles  gj  ’ computed  from  the  solution  f.  should  deviate  from  the  input 
data  profiles  g.  by  an  amount  equal  to  the  estimated  rms  errors,  e,  in 
g(z).  That  is 


No^  = [gj  - = Ne' 

j=l 


{B-6) 


The  smoothing  parameter  y in  (B-3)  is  adjusted  until  the  condition  (B-6) 
is  met.  The  smoothing  matrix  H also  incorporates  the  symmetry  and  end- 
point conditions 


on  the  smoothed  function. 


^0  ■ h 


^N+2  " ^Wl  ■ 


(B-7) 
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For  the  present  application,  the  method  must  be  modified  to  account 
for  two  additional  conditions.  First,  an  iterative  scheme  must  be  included 
in  order  to  account  for  nonthin  optical  depths,  and  second,  account  must 
be  made  of  the  fact  that  there  is  not  just  one  integral  equation  to  be 
inverted,  but,  in  fact,  two  coupled  equations.  The  two  equations  are 
the  radiance  and  absorptance  equations  [see  (81)  and  (82)  of  main  text)] 


N(z)  = 2 ^ J(r)  w‘^(z.r) 


2 2 
r - z 


(B-8) 


W(2)  = 2 ^ K(r)  y°(z,  r) 


2 2 
r - z 


(B-9) 


^(z)  and  W(2)  are  the  function  g(z)  in  (B-1)  and  J(r)  and  K(r)  are  the 
unknown  function  f(r),  2w®(z,  r)/(r^  - z^)^^^  and  2y^{z,r)/(r^  - z^)^*^^ 
are  the  kernel  function  K(2,  r)  of  (B-1),  and  a = z,  b = R.  In  the  solution 
of  (B-1),  K(z,  r)  is  assumed  to  be  known.  The  kernels  of  (B-8)  and 
(B-9),  however,  are  not  known  because  the  functions  w^(z,  r)  and  y^(z,  r) 
are  unknown.  As  for  the  iterative  Abel  inversion  method,  however,  we 
do  know  that  for  optically  thin  sources  w = y = 1 and  we  can  effect  an 
iteration  starting  with  these  values.  With  the  c(r)  and  T(r)  profiles 
obtained  with  these  initial  values  for  w^  and  y^,  new  values  for  w^  and  y^ 
can  be  computed,  and  the  process  continued  until  convergence  is  obtained. 
This  iteration  scheme  is  diagramed  in  Fig,  B-1. 

The  second  consideration  is  that  instead  of  a single  constraint  of  the 
form  (B-6),  there  are  two,  and  accordingly,  there  are  two  smoothing 
parameters  y that  must  be  varied  until  the  constraints  are  simultaneou.3ly 
met.  This  requirement  of  double  variation  is  a serious  defect  of  the 
method  since  a lot  of  computer  time  can  be  spent  searching  for  the  proper 
values  of  y. 

The  Phillips -Twomey  method  was  applied  to  the  moderately -deep 
inversion  example  considered  in  Section  VUB  and  VIIIB  of  the  main  text. 


-223- 


Fig.  B-1.  Phillips-Twomey  Inversion  Flow  Diagram 


When  applied  to  computed  transverse  E/A  profiles  with  no  superimposed 
random  errors,  the  convergence  chc^racteristics  are  very  similar  to 
those  of  the  iterative  Abel  inversion  method.  Comparisons  between  the 
results  of  the  iterative  Abel  and  Phillips-Twomey  inversion  methods  for 
the  case  of  five  percent  superimposed  random  errors  are  shown  in  Fig. 
B-2.  The  Phillips-Twomey  smoothing  parameters  were  taken  as 
■Yg,  = = EOOO  in  order  to  achieve  a values  near  five  percent.  Except 

for  values  near  r/R  = 1,  the  Phillips-Twomey  method  is  in  essential 
agreement  with  the  iterative  Abel  method.  The  Phillips-Twomey  results 
appear,  in  fact,  to  be  a somewhat  smoothed  version  of  the  iterative  Abel 
results. 
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APPENDIX  C 

WIDE  BAND  MODEL  PARAMETERS 


Within  the  random  band  model  formulation  for  an  array  of  Lorentz 
lines  in  a sx>ectral  interval  Av  centered  on  v,  the  mean  transmittance  in 
Av  for  a uniform  optical  path  is 


t(v) 


-P(v)F[x(v)] 


k(v) 

P(v) 


, . k(v)u 

fundamental  band  model  parameters  (known  for  each 
Av  over  a whole  vibration  rotation  band) 


(C-1) 


u = optical  depth 

F(x)  = curve -of -growth  function - 


V 2x/it 


Direct  use  of  this  expression  for  t(v)  in  an  experimental  application  pre- 
supposes detection  with  a rectangular  instrument  function  of  width  Av. 

Consider  detection  over  a band  pass  v — ► v,.  Let  the  external  radia- 

1 w 

tion  source  (with  which  the  absorption  measurements  are  made)  have  a 
spectral  output  S(v);  the  detector  (less  any  filter)  response  be  R(v);  and  the 
filter  transmission  be  T(v),  The  total  intensity  detected  is 


I = ^ S(v)T(v)R(v)T(v)dN 
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Without  the  intervening  gas,  the  intensity  detected  would  be 

''2 

Iq  = ^ S(v)T(v)R(v)dv 

The  apparent  wide  band  transmittance  of  the  gas  is  then 

^'2 

T = C o(vK(v)dv  (C-2) 

where 

, .._SCvJTiy)Riv),_. 

^2 

S(v)T(v)R(v)dv 


Substitution  of  (C-1)  into  (C-2)  gives  the  exact  wideband  transmittance 
function 


2 

Use  of  this  transmission  function  in  practical  applications  is  quite  incon- 
venient. A*’  an  approximation  to  (C-3),  assume  that  T can  be  expressed 
in  the  same  form  as  (C-1)  but  in  terms  cf  effective  parameters  and 

P . That  is,  assume 
e 
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- ^ ^-PeF(xe) 


k u 
e 


(C-4) 


The  k and  6 are  to  be  deternained  bv  forcing  (C-3)  and  (C-4)  be 
e e ' 

equal  in  sonae  limits  of  absorption.  Between  these  limits,  (C-4)  is 
assumed  to  be  a good  approximation  to  (C-3).  Setting  (C-3)  equal  to 
(C-4)  gives 


g-PeF(xe) 


"2 

^ Q'(v) 


p(v)F[x(v)]  , 
e dv 


T 


(C-5) 


First  consider  the  so-called  weak-  and  strong-line  absorption  limits. 
These  obtain,  respectively,  as  x — ►Oandx — ► ■».  For  x — ►O,  (C-5) 
reduces  to 


1 - p X = r a(v)[l  - p(v)x(v)]dv 


e e 


k u = ^ Q(v)k(v)udv 


and  thus  provides  the  definition  for  k 


An  important  point  must  be  made  clear  here  These  limits  refer  to  how 
each  line  in  the  array  absorb.  Even  though  each  line  is  absorbing  strongly, 
the  absorptance  effected  by  the  array  may  still  be  very  small  if  the  lines 
are  well  isolated  from  each  other  (compared  to  their  width).  Thus,  the 
strong  line  absorption  limit  does  not  necessarily  correspond  to  a strong 
band  absorption  limit. 
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(C-6) 


^ Q!(v)k(v)dv 


As  X — ►<»,  (C-5)  reduces  to 


which,  when  solved  for  p , yields 

e 


2uk 


-in 


e 


\ Qf(v)e 


Vfp' 


v)k(v)u 


dv 


{C-7) 


Equations  (C-6)  and  (C-7)  provide  definitions  for  the  effective  band 
model  parameters  to  be  used  in  (C-4)  when  the  criterion  for  matching 
(C-4)  to  {C-3)  is  correspondence  in  the  limits  of  weak-  and  strong-line 
absorption. 


Two  other  traditional  absorption  limits  are  for  overlapping  and  non- 
overlapping lines.  These  limits  are  obtained,  respectively,  for  p — ► <*>,  and 
p — ►O.  For  p — ►«,  X — ►O  and  (C -5)  reduces  to 


and  yields  the  definition  for  k 
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k 

e 


- In 
u 


For  P — ►O,  on  the  other  hand,  (C-5)  reduces  to 


(C-8) 


/ - k P u = C Q'(v) 

1 - \/— k(v)p(v)u 

/ IT  e e 

V ^ 

dv 


and  yields  the  definition  for  P 


r'^2. 


^ Qf(v)''  k(v)P(v)  dv 


(C-9) 


Which  two  of  the  four  absorption  limits  are  chosen  to  define  k and 
P^  should  be  dictated  by  the  experimental  conditions  of  the  application  of 
the  model.  Consider  the  Transtage  conditions.  For  a line -of-sight  through 
a full  plume  diameter,  the  approximate  conditions  of  Table  C-1  obtain 
for  the  NERD  and  (narrow)  bandpasses. 

The  first  observation  from  this  table  is  that  the  condition  of  nonover- 
lapping lines  prevails.  For  H O,  y/b  0.  0013  and  for  CO  , y/b  0.  043. 

£•  Ct 

This  result  strongly  suggests  the  use  of  (C-9)  for  the  definition  of  fi^. 

The  second  observation  is  that  x never  obtains  a value  so  large  that  the 
strong  line  absorption  limit  can  be  said  to  have  been  achieved.  The  table 
entries  are  for  the  full  line  of  sight  through  the  plume.  In  order  to  per- 
form radiance  calcxilations,  t must  be  knowTi  for  all  distances  into  the 
plume.  Although  the  strong  line  absorption  limit  does  not  obtain  for  the 
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Table  C-1.  Approximate  Transtage  Conditions 


L 120  cm 
T “•  1400  K 
p 0.  013  atm 


Variable 

^°2 

c 

0.4 

0.  04 

u = cpL  (atm  cm) 

0.  62 

0.  062 

k (cm  ^ / atm) 

0.  06 

3 

1/6  ( 1 / cm  ^ ) 

2 

100 

Y/p  (cm  Vatm) 

0.05 

0.033 

p = ^ 

6 

0.  008 

0.  27 

X 

4.  7 

0.  69 

F(x)  = ~|^v/  1 + 2itx  - 1 J 

1.4 

0.42 

PF(x) 

0.  012 

0.  11 

-PF(x) 

T = e 

0.99 

0.89 

.•!f^Si-.js*H*5c  . 


full  distance  through  the  plume,  the  weak  line  absorption  limit  must  obtain 

for  small  distances  into  the  plume.  Therefore,  (C  6)  is  an  appropriate 

definition  for  k . 

e 

For  the  Tranctage  conditions,  comparisons  between  (C-7)  auad  {C-9) 
for  gave  nearly  identical  results.  [Note,  (C-7)  was  applied  with  an  arbi- 
trarily large  u aind  found  to  give  0^.  nearly  independent  of  u,]  This  result  is 
easily  explained  by  noting  that  the  exponential  argument  of  (C-?)  is  only 
~0.0i4  for  H^O  auid  ~0.  18  for  CO2.  Then,  first  order  expansions  of  the 
exponential  and  logarithm  are  justified,  and  the  resulting  expression  is  (C-9). 
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Facilities 


AFRPL 


AEDC 


NASA 


ACRONYMS 


- Air  Force  Rocket  Propulsion  Laboratory 

- Arnold  Engineering  and  Development  Center 

- National  Aeronautics  and  Space  Administration 


AFGL  ~ Air  Force  Geophysics  Laboi'atory 


Test  Programs 


HAPRAP  - High  Altitude  Plume  Radiation  Program 
ERASE  - Emitted  Radiation  from  Special  Engines 


NERD 


- Nozzle  Exit  Radiation  Determination 


Technir-'.l 


E/A  - emission  and  absorption 

pTc  - pressure,  temperature  and  species  concentration 

EMABIC  - Emission/Absorption  inversion  Code 

CG  - Curtis-Godson  nonuniformity  approximation 

DR  - derivative  nonuniformity  approximation 

LS  - Lindquist- Simmons  nonuniformity  approximation 

ML  - mean  line  nonuniformity  approximation 

RW  - Rodgers  and  Williams  Voigt  model 


FOV 


field  of  view 


- oxidizer -to -fuel  ratio 


LABORATORY  OPERATIONS 

Th«  Lftboratory  OperAtions  of  The  Aerotp&ce  Corporation  it  conducting 
Mperimentnl  end  thcoreticel  inveetigetioD*  necopnery  for  the  eveiue&ion  end 
eppUcetlon  of  ecientiflc  edvencee  to  new  miUtery  concepte  end  eyetene.  Ver* 
•etlUty  end  flexibility  heve  been  developed  to  e high  degree  by  the  Ihboretory 
pereonnel  in  deeling  with  the  meny  problems  encountered  in  the  netion'e  repidly 
developing  epece  end  miee^ie  eyeteme.  Exnertiee  in  the  leteet  ecientific  deveU 
opmente  le  yitel  to  the  eccon^liebznent  of  teeke  releted  to  theee  probleme.  The 
leboretoriee  thet  contribute  to  thie  reeeercK  ere: 


Aerophvice  L^iboretory:  lAuncb  end  reentry  eerodynemice^  heet  trece- 
fer«  ree^^  pi^yeCce.  cKem'cel  kinetice*  str^icturel  mecbenice,  flight  dynemice* 
etxnoepheric  pollutlotu  end  high-power  gee  leeere. 

Chemletry  end  Phyeice  Leboretory;  Atmospheric  reectione  end  etmcso 
pheric  ^tics,  ^emicei  reel^one  Tn  poUuted  etmosph  res,  cheitdcel  reactions 
of  excitt^  species  in  rocket  pluxres,  chemical  thermc<^ynemics,  pteeme  end 
leser-induced  reectione,  leeer  ch'^mlstry,  propulsion  chemletry,  spece  vecuum 
end  redietion  effects  on  meteriels.  lubricetion  end  suri'ece  phenomene,  photo- 
sensitive meteriels  end  sensors,  high  precision  leser  ranging,  end  the  appll- 
cetlon  of  physics  end  chemietr**  to  probleme  of  lew  ei^'^icement  end  biomedicine, 

Electronics  Research  Leboretory:  Electromagnetic  theory,  devices,  end 
propagation  pWnornen^  Including  pleeme  electromagnetics;  quantum  electronics, 
lasers,  end  electro-optics;  communication  sciences,  applied  electronics,  semi* 
conducting,  supercon^cting,  end  crystal  device  physics,  optical  end  ecoustlcel 
imaging;  atmospheric  pollution;  mllUmeter  wave  end  far-infrared  technology. 

Meteriels  Sciences  LAboretory:  Development  of  new  meteriels;  metal 
matrix  composites  en^new^  lor^ms  <n  carbon;  test  and  eveluetioB  of  graphite 
end  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  envlror,ment;  application  of  fracture  mechanics  to  stress  cor* 
roalon  and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion fr<^  t&e  atmosphere,  <Iensity  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetosphen-  ,)hysics,  coerce  rays,  generation  and  pri^agation 
of  plasma  waves  In  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
Reids:  space  astronomy,  x-tay  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solsr  activity  on  the  earth's  atmosphere,  ionosphere,  and 
m^netosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions ic  space  on  space  systems. 

THE  AEROSPACE  CORPORATION 
El  cegundo,  California 


